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SELF-GUIDING  AND  STABILITY  OF  INTENSE 
OPTICAL  BEAMS  IN  GASES  UNDERGOING  IONIZATION 

I.  Introduction 

The  propagation  of  optical  pulses  in  gases  is  relevant  to  a  wide  range 

1  2 

of  applications,  such  as  ultra-broadband  optical  generators,  ’  optical 
harmonic  generators, x-ray  lasers,^  and  laser-driven  accelerators 
For  these  applications  it  is  necessary  that  the  optical  pulse  be  intense 
and  propagate  extended  distances.  In  the  absence  of  an  optical  guiding 
mechanism  the  propagation  distance  is  limited  to  approximately  a  Rayleigh 
(diffraction)  length.  At  sufficiently  high  power  and  intensity  the 
propagation  distance  is  strongly  affected  by  nonlinear  self-focusing  and 
ionization  (plasma  generation). 

An  optical  beam  propagating  in  a  neutral  gas  is  affected  by 
diffraction,  refraction,  nonlinear  self-focusing,  ionization,  and  plasma 
defocusing.  Self-focusing,  for  example,  is  due  to  the  intensity  dependent 
part  of  the  refractive  index  and  occurs  when  the  optical  power  is  above  the 

•IQ  <1  ^ 

nonlinear  focusing  power.  As  the  beam  focuses,  the  increased 

intensity  results  in  ionization  and  plasma  formation  which  tends  to  defocus 

the  optical  beam^’^^”^^,  see  Fig.  1.  A  balance  between  the  nonlinear 

focusing  and  plasma  defocusing  can  result  in  a  self-guided  optical  beam. 

In  this  paper  the  propagation,  self-guiding,  and  stability  of  two 

types  of  optical  beams  are  analyzed.  The  two  beams  considered  are  a 

2  2 

fundamental  Gaussian  beam  of  the  form  E^exp(-r  /r^  +  it(»)e^/2  +  c.c.  and  a 

2  2 

higher-order  radially— polarized  beam  of  the  form  E^('r2r/r^)exp(-r  /r^  + 

i\|»)ej./2  +  C.C.,  where  E^  is  the  electric  field  amplitude,  r^  is  the  spot 

size  and  is  the  phase.  The  results  of  this  paper  include  (1)  envelope 

equations  describing  the  evolution  of  the  optical  beam  spot  size,  which  are 

22  23 

derived  by  using  the  source-dependent  expansion  method,  ’  (2)  the 

critical  power  for  nonlinear  self-focusing  of  the  higher-order  mode,  which 
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is  four  times  greater  than  that  of  the  fundamental  Gaussian,  (3)  self- 
guided  beam  solutions,  which  result  from  a  balance  of  nonlinear  self- 
focusing  and  plasma  defocusing,  (4)  the  analysis  of  a  new  ionization- 
modulation  instability,  which  disrupts  self-guided  beams,  (5)  the  evolution 
of  the  optical  beam  phase  velocity,  which  is  less  than  the  speed  of  light 
for  a  self-guided  beam,  and  (6)  a  new  configuration  of  an  inverse  Cherenkov 
accelerator,  which  is  based  on  a  self-guided,  radially-polarized,  higher- 
order  Gaussian  beam. 

One  important  application  of  intense  optical  pulses  propagating  in 
gases  is  laser-driven  electron  accelerators,  which  are  referred  to  as 
inverse  Cherenkov  accelerators  (ICAs).^”^  In  the  conventional  ICA,^  ^  the 
optical  beam  driver  can  consist  of  either  a  radially-polarized  i)  higher- 
order  Gaussian  mode  or  ii)  a  nonideal  first-order  Bessel  mode.  Associated 
with  these  modes  is  an  accelerating  axial  field  peaked  along  the 
propagation  direction  (z-axis).  In  general,  the  electron  acceleration 
distance  is  limited  by  either  the  diffraction  distance  or  the  electron 
slippage  distance.  Since  the  optical  beam  in  the  ICA  propagates  in  gas, 
the  phase  velocity  can  be  less  than  the  speed  of  light  and  controlled  by 
varying  the  gas  density.  Electron  slippage  is  minimized  by  matching  the 
electron  velocity  to  the  phase  velocity  of  the  accelerating  field.  The 
acceleration  distance,  however,  is  still  limited  by  the  diffraction  length. 
For  a  higher-order  Gaussian  mode  the  diffraction  length  is  a  Rayleigh 
length  which  is  precisely  the  slippage  distance  in  vacuum.  In  this  case, 
as  far  as  energy  gain  is  concerned,  there  is  essentially  no  advantage  in 
introducing  a  gas  since  the  effective  acceleration  length  is  limited  to 
approximately  a  Rayleigh  length  as  discussed  in  Sec.  V.  For  a  fixed  total 
optical  beam  power,  however,  the  energy  gain  in  the  conventional  ICA  can  be 
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significantly  increased  by  using  a  nonideal  Bessel  beam,  as  discussed  in 

Appendix  D.  To  further  enhance  the  energy  gain,  the  ICA  requires  self- 

9  6—8 

guiding  of  the  optical  driver.  Previous  studies  of  the  ICA  also 
neglect  the  intensity  dependent  effects  in  the  refractive  index,  i.e., 
nonlinear  self-focusing,  as  well  as  ionization,  i.e.,  plasma  effects.  We 
propose  and  analyze  a  self-guided  ICA  configuration  that  operates  at  laser 
powers  near  the  nonlinear  self-focusing  power  and  at  intensities  high 
enough  to  slightly  ionize  the  gas. 

Another  possible  application  of  intense  optical  beams  in  gases  is  the 

12  3  4 

generation  of  ultra-broadband  ’  or  harmonic  radiation.  ’  A  short  pulse 

optical  beam  propagating  in  a  nonlinear  medium  will,  among  other  things, 

13—16 

undergo  self-phase  modulation  which  results  in  frequency  broadening. 

Since  the  degree  of  frequency  broadening  increases  with  both  propagation 

distance  and  optical  intensity,  the  self-guiding  of  a  short  optical  beam 

may  be  well  suited  for  ultra-broadband  radiation  generation.  A  self-guided 

3  4 

Gaussian  beam  may  have  application  to  harmonic  generation,  ’  since  the 
propagation  medium  consists  mainly  of  a  neutral  gas  and  a  very  narrow 
plasma  column  along  the  axis.  The  harmonics  could  be  guided  by  the  driving 
optical  beam  and  phase  patching  may  be  achieved  by  introducing  a  background 
plasma. 

This  paper  is  organized  as  follows.  The  propagation  model  is 
presented  in  Sec.  II,  and  includes  discussions  of  the  wave  equation  in  a 
gas  undergoing  ionization,  the  linear  and  nonlinear  polarization,  plasma 
generation,  the  reduced  wave  equation,  photo-ionization,  the  solution  of 
the  wave  equation  using  the  source  dependent  expansion  method,  and  the 
resulting  equations  describing  the  evolution  of  the  envelope,  amplitude, 
and  phase  of  the  optical  beam.  The  propagation  of  a  fundamental  Gaussian 


3 


beam  is  examined  analytically  and  numerically  in  Sec.  Ill,  including  the 
envelope  equation,  self-guided  solutions,  and  the  stability,  i.e.,  the 
ionization-modulation  instability.  Numerical  results  on  the  propagation  of 
the  fundamental  Gaussian  beam  are  also  presented  in  Sec.  III.  The 
propagation  of  a  higher-order  radially-polarized  beam  is  analyzed  in  Sec. 
IV,  including  the  envelope  equation,  self-guided  solutions,  and  the 
stability.  Numerical  results  on  the  propagation  of  the  higher-order 
radially-polarized  beam  are  also  presented  in  Sec.  IV.  The  analysis  of  a 
self-guided  inverse  Cherenkov  accelerator  is  presented  in  Sec.  V. 
Attenuation  of  the  optical  beam  due  to  electron  collisions  and  ionization 
losses  is  analyzed  in  Sec.  VI.  Section  VII  contains  a  discussion  and 
summary.  This  paper  also  includes  three  Appendices  discussing  (A)  photo¬ 
ionization  rates,  (B)  electron  collision  frequencies,  (C)  radiative  and 
collision  losses  on  accelerated  electrons,  and  (D)  inverse  Cherenkov 
acceleration  with  Bessel  (axicon)  beams. 

II.  Propagation  Model 

The  propagation  of  intense  optical  beams  in  gases  is  affected  by  a 

combination  of  diffraction,  refraction,  and  ionization.  The  refractive 

13—16 

index  of  a  gas  generally  has  an  intensity  dependent  part,  n  =  n^  + 

n-I,  where  n  is  the  linear  refractive  index,  n«  is  the  nonlinear 

refractive  index  and  I  is  the  intensity  of  the  optical  beam.  Generally,  TI2 

is  positive  and  results  in  self-focusing  of  the  optical  beam  if  the  power 

is  greater  than  the  nonlinear  focusing  power.  The  nonlinear  focusing 

13—16  2 

(critical)  power  for  a  fundamental  Gaussian  beam  is  =  X  /(2mi^n2), 

where  X  is  the  vacuum  wavelength.  As  the  beam  self-focuses  the  peak 
intensity  increases  resulting  in  ionization  and  the  generation  of  a  plasma. 
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In  the  region  of  the  plasma  the  refractive  index  is  n(r)  =  n^  +  n2l  - 
92  2  1/2 

<0  (r)/2w  ,  where  »  =  (4iiq  n  /m)  is  the  plasma  frequency,  n^  is  the 
P  P  P  P 

plasma  density,  and  »  =  2llc/X  is  the  frequency  of  the  optical  beam.  The 

local  decrease  in  the  refractive  index  due  to  the  plasma  tends  to  defocus 
Q  17_21 

the  optical  beam.  ’  If  diffraction,  self-focusing  due  to  n2  and 

defocusing  due  to  plasma  generation  are  properly  balanced,  a  self-guided 

optical  beam  can  be  formed  and  propagated  over  extended  distances,  i.e., 

9  20  21 

many  vacuum  Rayleigh  lengths.  *  ’ 

Our  propagation  model  includes  a  number  of  assumptions.  The  short 
pulse  optical  beam  is  assumed  to  be  adequately  described  by  a  single  source 
dependent  Laguer re-Gauss ian  mode,  which  is  a  superposition  of  many  vacuum 
Laguerre-Gaussian  modes.  The  model  is  not  valid  when  the  optical  power 
greatly  exceeds  the  nonlinear  focusing  power,  since  the  beam  is  expected  to 
filament  into  higher-order  modes.  Ionization  is  considered  in  the  high 
field  limit^  ~  (Keldysh  parameter  less  than  unity)  and  is  modeled  by  the 
tunneling  ionization  rate,  see  Appendix  A.  The  attenuation  of  the  optical 
beam  due  to  ionization  and  collisional  losses  is  estimated  and  found  to  be 
small  enough  to  neglect.  The  nonlinear  polarization  field  of  the  gas  is 
included  to  third  order. in  the  optical  field  whereas  the  plasma  current  is 
included  to  first  order. 

A.  Wave  Equation  in  Gas  Undergoing  Ionization 

The  dynamics  of  optical  beams  propagating  in  a  gas  undergoing 
ionization  is  governed  by  the  wave  equation, 

-  c"^3^/3t^)E  =  4nc”^(3^P/3t^  +  3Jp/3t),  (1) 

9  2  2  2. 

where  E  is  the  electric  field  of  the  optical  beam,  v  =  +  3  /3z  ,  z  is 

the  axial  propagation  direction,  P  is  the  polarization  field  associated 
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with  the  gas  and  Jp  is  the  plasma  current  density  associated  with  the 
ionized  gas.  In  obtaining  Eq,  (1)  we  have  neglected  a  small  source  term 
proportional  to  the  gradient  of  the  plasma  density. 


1.  Linear  and  Nonlinear  Polarization 

The  polarization  field  can  arise  from  a  number  of  processes;  these 

include  electronic  polarization,  molecular  orientation,  electrostriction, 

saturated  absorption  and  thermal  effects.  In  the  present  paper  we 

will  be  concerned  with  changing  the  refractive  index  on  a  fast  time  scale, 

-12 

typically  less  than  10  sec.  On  this  time  scale  the  electronic 
polarization  is  the  dominant  contribution  to  the  nonlinear  refractive  index 
and  is  due  to  the  optical  field  modifying  the  atomic  electronic 
distribution. 

In  the  simple  Lorentz  model  of  the  atom  the  electrons  are  assumed 

to  consist  of  a  charge  distribution  oscillating  in  an  effective  potential. 
Nonlinearities  in  the  effective  potential  result  in  a  field  dependent 
refractive  index  for  the  medium.  In  the  following  description  of  the 
polarization  field  only  isotropic  matter  having  ensemble  averaged  inversion 
symmetry,  i.e.,  cen t rosy mme trie  ensemble  averaged  effective  potentials, 
will  be  considered.  This  includes  all  liquids,  gases,  amorphous  solids  as 
well  as  many  crystals. 

The  electric  polarization  field  is  defined  by  P  =  qn^^,  where  q  is 
the  electronic  charge,  n^  is  the  density  of  atoms  or  molecules  and  x^  is 
the  displacement  of  the  electronic  distribution  from  equilibrium  due  to  the 

optical  field.  The  polarization  field  in  the  classical  single  resonant 

.  j  ■  •  K  13-16 

frequency  model  is  given  by 

2r3P/3t  =  (q^n^/m)E,  (2) 
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where  is  the  characteristic  resonant  frequency  of  the  electronic 
distribution,  is  a  constant  associated  with  the  nonlinear,  i.e., 
nonparabolic,  nature  of  the  effective  potential,  is  a  normalizing 
polarization  field  amplitude  and  P  is  a  damping  term.  Equation  (2)  is  an 
accurate  description  for  the  polarization  field  when  the  optical  frequency 
is  far  from  the  resonant  frequencies.  Typically,  the  resonant  frequency  ^ 
is  in  the  ultra  violet  regime,  »  w.  The  polarization  field  given  by 
Eq.  (2)  contains  dispersion,  damping  and  third  order  nonlinear  effects. 

In  the  limit  where  i)  dispersive  effects  are  weak  (far  from  resonance. 


»  <o) ,  ii)  damping  effects  can  be  neglected  (F  <<  «),  and  iii)  the 
nonlinear  term  in  Eq.  < 
can  be  approximated  by 


nonlinear  term  in  Eq.  (2)  is  small  (S^P^/P^  «  S2^),  the  polarization  field 


P  =  +  X^^^<E-E>E 


=  (l/4n)(n^  -  1  +  2n^n2l)E,  (3) 

where  X^^^  =  q^n^/(inS^)  «  1  is  the  constant  linear  susceptibility,  X^^^  = 

«  X^^^/<E-E>  is  the  constant  third  order  susceptibility 
of  the  neutral  gas,  the  brackets  <  >  denote  a  time  average,  n^  =  (1  + 
is  the  linear  refractive  index  of  the  neutral  gas,  ^^2  ~ 
(8ii^/n^c)X^^^  is  the  nonlinear  component  of  the  refractive  index,  I  = 
(c/4n)n^<E-E>  is  the  intensity  and  |n2l  |  «  has  been  assumed. 


2.  Plasma  Generation 

The  ionization  of  the  background  gas  by  the  optical  beam  results  in 

the  generation  of  plasma  electrons.  The  plasma  current  density  is  given  by 

J  =  an  V  ,  where  n  and  v  are  the  plasma  density  and  fluid  velocity 
~P  P~P  P  ~P 
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respectively.  To  lowest  order  in  v^,  the  continuity  and  fluid  velocity 
equations  are 


3np/3t  +  ^‘(OpVp)  =  S, 

(4a) 

mn  3v  /3t  =  qn  E  -  mv  S, 

p  ~p  p~  ~p 

(4b) 

where  S  is  the  plasma  source  term  proportional  to  the  ionization  rate,  E  is 

the  optical  electric  field,  the  v^  x  B  force  and  thermal  effects  are 

neglected  in  Eq.  (4b),  and  the  electrons  are  assumed  to  be  created  with 

zero  velocity  when  ionized.  Combining  Eqs.  (4a, b),  and  keeping  terms  to 

.  .  .  .  U  29-31 

lowest  order  in  v^,  the  plasma  current  density  is  given  by 

3Jp/9t  =  (Wp/4n)E,  (5) 

where  w  =  (4nq^n  /m)^^^  is  the  electron  plasma  frequency.  The  evolution 
P  P 

of  the  plasma  density  depends  on  the  photo-ionization  rate  and  is  discussed 

later.  In  obtaining  Eq.  (5)  nonlinear  and  collisional  effects,  see 

Appendix  B,  have  been  neglected.  Ionization  and  collisional  losses  are 

analyzed  in  Sec.  VI  and  found  to  be  small.  Nonlinear  plasma  effects  are 

small  compared  to  the  nonlinear  neutral  gas  effects  which  are  represented 

by  the  term  n2l  in  Eq.’  (3).  The  magnitude  of  nonlinear  plasma  effects 

compared  to  nonlinear  neutral  gas  effects  is  approximately  given  by  the 

11  32 

ratio  of  the  critical  power  for  relativistic  focusing  ’  to  the  nonlinear 
focusing  power  and  is  found  to  be  negligibly  small. 


3.  Reduced  Wave  Equation 

The  propagation  of  the  optical  beam  is  described  by  Eqs.  (1),  (3)  and 
(5)  together  with  the  tunneling  ionization  model  discussed  in  Sec.  II. B. 

To  proceed  with  the  analysis,  it  is  convenient  to  transform  from  the  (z,t) 
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coordinates  to  the  (t, ti)  coordinates,  where  ^  =  z  -  vt  and  h  =  z*  For  a 
beam  propagating  in  the  positive  z  direction  with  group  velocity  v,  ^  is 
the  distance  behind  the  front  (K  =  0)  of  the  optical  beam  and  h  is  the 
propagation  distance.  In  these  new  coordinates,  the  optical  field  has  the 
form 

E  =  E  exp(ik,t)/2  +  c.c.,  (6) 

where  E(r, 0,^,1!)  is  the  complex  amplitude  and  is  a  slowly  varying  function 
of  K  and  k  =  «/v,  co  =  2iic/X  is  the  optical  frequency,  X  is  the  vacuum 
wavelength,  and  c.c.  denotes  the  complex  conjugate.  In  the  (^,Vi) 
coordinates,  Eq.  (1)  can  be  written  in  the  paraxial  approximation  as 

(vj  +  2ik3/an)E  =  K^E,  (7) 

2 

where  K  is  given  by 

-  2k^(n.,/n  )I.  (8) 

P  z  o 

In  obtaining  Eq.  (7)  we  used  the  transformations  3/3z  =  3/^  +  3/3>i  and  3/3t 

=  -  v3/3^.  The  linear  group  velocity  as  well  as  the  linear  phase  velocity 

2 

is  V  =  w/k  =  "Fhe  paraxial  approximation  implies  that  K  is  small 

2 

compared  to  k  .  In  the  absence  of  the  nonlinear  index  (n2  =  0),  the  , 

paraxial  approximation  requires  that  the  plasma  density  be  small  compared 

2  2 

to  the  critical  density,  Wp  <  «  . 

B.  Photo-Ionization  Model 

33-35 

Ionization  can  occur  by  electron  collisional  processes  or  by  the 

0  y  0  Q 

intense  optical  fields  directly,'^  ”  i.e.,  photo-ionization.  In  the 

absence  of  collisions  or  for  laser  pulses  short  compared  to  a  collision 
time,  photo-ionization  is  the  dominant  process.  Photo-ionization  can  take 
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place  by  either  tunneling  or  multi-photon  processes,  see  Appendix  A.  These 

1/2 

regimes  are  characterized  by  the  Keldysh  parameter  Yj^  =  ’  where 

2 

U,  is  the  ionization  energy  and  e  =  (l/2)m(q |E |/m«)  is  the  electron 

X.  OS 

oscillation  energy.  The  low  field  limit  (yj^.  >  1)  corresponds  to  the  multi¬ 
photon  ionization  regime,  whereas  the  high  field  limit  (yj^  ^  1)  corresponds 
to  the  tunneling  ionization  regime. 

The  evolution  of  plasma  density  in  Eq.  (4a)  is  given  by 


9np/8t  =  -  np)W(|i|),  (9) 

where  n  is  the  initial  neutral  density,  W(|e|)  is  the  ionization  rate  and 
no 

the  convection  term  y*(n  v  )  is  neglected.  For  a  linearly  polarized  optical 

p-p 

24-28 

field,  the  ionization  rate  in  the  tunneling  limit  (yj^  <  1)  is  given  by 

W(|El)  =  4(3/n)^^2Q^(Uj)^^‘^(Eg/|E|)^^2exp(-  |  |E  |]  , 

(10) 


16  1 

where  fi  =  a,c/a  =  4.1  x  10  sec~  is  the  characteristic  atomic 
0  r  o 

-9 

frequency,  =  1/137  is  the  fine  structure  constant,  a^  =  5.3  x  10  cm  is 
the  Bohr  radius,  Uj  =  U^/Ug,  is  the  ionization  energy  in  eV,  Ujj  =  13.6 
eV  is  the  ionization  energy  of  hydrogen,  |e|  is  the  magnitude  of  the 
optical  field  and  Eg  =  Iqj/a^  =  5.2  GV/cm  is  the  atomic  field  of  hydrogen. 
The  intensity  of  a  linearly  polarized  Gaussian  optical  beam  in  vacuum,  with 
a  peak  field  equal  to  Eg,  is  I  =  (c/8il)Eg  =  3.6  x  10^^  V/cm^.  Equation  (9) 
assumes  that  the  gas  is  at  most  singly  ionized.  The  solution  of  Eq.  (9) 
yields 


k 


2 

P 


exp 


o 


[-(n^/c)  J  W(|E|)d^') 


K 


(11) 
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where  we  have  set  8/3t  =  -(c/n^)9/3£;  in  Eq.  (9),  ^  is  defined  in  the  region 

^  <  0,  i,  =  0  corresponds  to  the  front  of  the  beam  and  ck  =  «  = 

—  po  po 

2  1/2 

(4nq  n^^/m)  is  the  plasma  frequency  associated  with  the  initial  neutral 

gas  density.  For  low  levels  of  ionization,  i.e.,  n^  «  n^^,  Eq.  (11) 
reduces  to 

o 

''p  ' 

t 

The  weakly  ionized  limit  is  sufficient  to  describe  self-guiding  of  optical 
beams,  since  it  will  be  shown  that  in  the  highly  ionized  limit  n^  =  n^^ 
there  is  no  matched  beam  solution.  The  expression  in  Eq.  (8)  can  be 
written  as 

o 

=  kpo(n^j/c)  J  W(|E|)d^'  -  2k2(n2/n^)I.  (13) 

C.  Source  Dependent  Expansion  Method 

The  following  analysis  is  based  on  the  SDE  method  developed  in  Ref. 

22.  The  SDE  is  a  powerful  method  for  solving  the  paraxial  wave  equation 
for  optical  beams  propagating  in  nonlinear  media.  In  the  SDE  method,  the 
optical  beam  is  expanded  in  a  complete  set  of  source  dependent  orthogonal 
Laguerre-Gaussian  functions.  These  functions  are  implicitly  functions  of 
the  propagation  distance,  >1>  through  the  optical  beam  parameters,  i.e., 
spot  size,  wavefront  curvature,  amplitude  and  phase.  The  optical  beam  can 
be  described  by  four  coupled  first  order  differential  equations  for  the 
beam  parameters  in  the  variable  h-  In  general,  E  can  be  written  in  terms 
of  a  complete  set  of  Laguerre-Gaussian  functions,  i.e.,  source-dependent 
modes. 
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(14) 


1=  X  ^ 

m,p 

where  m,p  =  0,1,2,...,  a„  ^(n),  b„  ^(H)  are  complex  coefficients  and  are 

m,  p  ,  p 

functions  of  ti>  ej  denote  transverse  unit  vectors  defining  the 
polarization,  dP(X)  *  xP'^^LP(X)exp(-  (1  -  ia)X/2),  X  =  Zr^/r^,  r^Ch)  is 
real  and  denotes  the  spot  size,  a(t|)  =  kr  /(2R  )  is  real,  R  is  the  radius 

SC  c 

of  curvature  associated  with  the  wave  front,  and  lP  is  an  associated 
Laguerre  polynomial,  e.g.,  Lq(X)  =  1  and  L^(X)  =  1  +  p  -  X*  The 
representation  in  Eq.  (14)  forms  a  complete  set  and  can  be  used  to 
represent  any  arbitrary  optical  beam. 

To  proceed  with  the  SDE  analysis  we  substitute  Eq.  (14)  into  the 
paraxial  wave  equation,  Eq.  (7),  carry  out  the  indicated  differential 
operations,  perform  the  operation 


2li 

J  ^cos(p'6),  sin(p'0) jd0/2ii 
o 


n  ^ 

on  both  sides  of  the  equation,  multiply  both  sides  by  (dP(X))  and  finally 
integrate  over  X  from  0  to  ®.  The  algebraic  details  can  be  found  in  Ref. 
22.  The  resulting  equation  for  a  is 

m,  p 


+  A 


m,p;  m,p 


imBa  .  -i(m  +  p  + 
m-l,p 


l)B*a. 


m+l,p 


-  iF  , 

m,p’ 


(15) 


where 

A  =  r  /r  +  l(2m  +  p  +  1) f(l  +  a^)/(kr^)  -  or  /r  +  0/2],  (16a) 
ni,pss  \  ss>£>/ 

B(h)  =  -  -  (1  -  a^)/(krg)  +  i/2  -  i  (r^/r^  -  2a/(krJ)] , (16b) 
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C30 


2n 

fm.p  -  2k  (i-TT)!  J  ■*«  J  ■'"(X.e.wicx.e.wej^ 

o  o 


X  (DP(X))*cos(pe)/(l  +  &p^o)» 


(16c) 


6  ,  is  the  Kronecker  delta,  the  dot  denotes  the  operator  3/9h,  and  the 

P»P 

asterisk  denotes  the  complex  conjugate.  The  equation  for  b  is  identical 

in,  p 

to  the  one  for  a  except  in  the  expression  for  F  ,  cos(pe)/(l  +  8  ) 

m,p  m,p  p,o 

is  replaced  with  sin(p0).  In  obtaining  Eq.  (15)  a  number  of  identities 
associated  with  Laguerre-Gaussian  functions  were  used,  including  the 
orthogonality  relation 


J  DP(X)(DP(X))*dX 

o 


n! 


m,n 


Equation  (15),  together  with  the  definitions  in  Eqs.  (16),  describes 
the  evolution  of  the  various  source  dependent  Laguere-Gaussian  modes. 
However,  Eq.  (15)  is  underdetermined  since  there  are  more  unknowns  than 
equations.  An  additional  constraint,  i.e.,  a  specification  of  the  function 
B(ti),  is  necessary  to  solve  Eq.  (15).  The  individual  source  dependent 

modes  in  Eq.  (14)  are  functions  of  the  spot  size,  r  (n),  wavefront  radius 

7. 

of  curvature,  =  krg/(2a),  amplitude  and  phase,  p.  Since  B(h)  is  also 

a  function  of  r  and  a,  the  evolution  of  the  source  dependent  mode  is 
s 

governed  by  the  particular  choice  for  the  function  B.  For  example,  if  we 
choose  B(h)  =  0,  we  recover  the  conventional  vacuum  modes.  In  general, 
expansion  in  terms  of  the  vacuum  modes  (B  =  0)  requires  many  modes  to 
accurately  describe  a  guided  optical  beam  over  distances  of  many  Rayleigh 
lengths.  A  more  appropriate  choice  for  B(q)  will  depend  on  the  particular 
problem  under  consideration  and  will  be  discussed  later. 


13 


The  dynamics  of  i)  a  fundamental  Gaussian  beam  and  ii)  a  higher-order 
radially-polarized  axially  symmetric  beam  will  be  considered.  The 
fundamental  Gaussian  beam  is  described  by  the  mode  numbers  m  =  0  and  p  =  0, 
whereas  the  higher-order  radially-polarized  axially  symmetric  beam  is 
described  by  the  mode  numbers  m  =  0  and  p  =  1.  The  analysis  can  be 
significantly  simplified  by  setting  m  =  0. 

In  the  following,  it  is  assumed  that  the  dynamics  of  the  optical  beam 
can  be  adequately  described  by  the  behavior  of  a  single  source  dependent 
mode,  e.g.,  the  m  =  0,  p  =  0  mode  for  the  fundamental  Gaussian  and  the  m  = 
0,  p  =  1  mode  for  the  higher-order  beam.  In  the  SDE  method,  it  is  assumed 
that  the  coupling  to,  as  well  as  the  amplitude  of,  the  higher-order  source 
dependent  modes  are  small.  In  fact,  an  optimal  choice  for  the  function 
B(n)  can  be  determined  from  Eq.  (15)  by  requiring  that  the  higher-order 
source  dependent  modes,  i.e.,  m  >  1,  are  small.  Since,  for  the  cases  under 
consideration,  |aQ  p|  »  la^j^  p|  for  m  >  1,  it  is  clear  from  Eq.  (15)  (with 
m  =  1)  that  the  optimal  choice  for  B  is 


B  = 


F  /  3.^  > 

l,p  0,p' 


(17) 


where  B  and  are  given  by  Eq.  (16b)  and  (16c)  respectively.  With  this 

1  >  P  . 

choice  for  B,  Eq.  (15)  (with  m  =  0)  yields 


^0,p)' 


■0,p 


(18) 


where  and  are  given  by  Eqs.  (16).  Equations  (17)  and  (18) 

0,p  U,p 

completely  determine  the  evolution  of  the  source  dependent  optical  beam 

mode.  Substituting  Eqs.  (16a, b)  into  Eqs.  (17)  and  setting  aQ^p  = 

E  exDfi©  ).  where  E  and  0  are  real  and  denote  the  field  amplitude  and 
o  o'  ^  o  o 

22 

phase,  we  obtain 
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(19a) 


+  (1  +  p)[(l  +  a^)/(kr^)  -  oi-g/r^  +  ^t/2]  =  -  G, 


or  /r  +  (1  -  a^)/(kr^)  -  i/2  =  -  H, 
s  s  s 


Tg/ts  -  2a/(krg)  =  0, 


(19b) 

(19c) 


E  /E  +  r  /r  =0, 
o  o  s  s 


(19d) 


where  G  =  /a^  and  H  =  F.  /a^  are  real.  The  source  functions  G  and 

0,p  0,p  If p  OfP 

H  are  given  by 


y  =  J  dx  K^(X,h)xPexp(-X) 


'  1 

a-x/(i+p)J 


(20) 


2 

and  K  is  given  by  Eq.  (13). 


D.  Envelope,  Amplitude  and  Phase  of  Optical  Beam 

Equations  (19b)  and  (19c)  can  be  combined  to  form  an  envelope  equation 
for  the  optical'  beam 

a^r^/aq^  -  4[l  +  krgH]/[k^rJ]  =  0.  (21) 

In  addition,  the  amplitude,  phase,  curvature,  and  axial  phase  velocity  v^j^ 
of  the  optical  beam  are  given  respectively  by 


a(E^rg)/aq  =  0, 

(22a) 

=  -d  +  p)  [2/(kr2)  +  h]  -  G, 

(22b) 

a  =  krg/(2R^)  = 

(22c) 

V  =  <1  - 

(22d) 
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where  |e/k|  «  1  was  assumed  in  Eq.  (22d).  Note  that  Eq.  (22a)  implies 

2 

that  the  optical  power,  which  is  proportional  to  ,  is  a  conserved 

2 

quantity,  consistent  with  the  paraxial  wave  equation  when  K  is  real. 

For  propagation  in  vacuum  (n^  =  1)  the  solution  of  Eqs.  (19a-d)  yield 

1  C  O  ^ 

the  conventional  vacuum  modes.  ’  In  vacuum,  the  source  terms  vanish, 

i.e.,  G  =  H  =  0,  and  the  solutions  are  characterized  by  a  spot  size  r^  = 

991 /9  * 

r  (1  +  n  /Z„  )  ,  r  is  the  minimum  spot  size  at  the  focal  point  q  =  z  = 

so'  Ro  so 

0,  =  krgQ/2  is  the  Rayleigh  length,  a  =  q/Zj^^  =  z/Zr^.  a  wavefront 

2  2 

radius  of  curvature  R^  =  z(l  +  Z^^/z  ),  a  phase  factor  6^  =  6^(11  =  0)  -  (1 
+  p)tan~^a,  an  amplitude  =  E^(q  =  0)r  /r  ,  and  a  phase  velocity  v  . /c  = 

O  O  so  5 

1  +  2(1  +  p)/k^rg. 

III.  Fundamental  Gaussian  Beam  Propagation 

We  first  consider  the  dynamics  of  a  fundamental  Gaussian  optical  beam 
propagating  in  a  gas  undergoing  ionization.  The  fundamental  Gaussian  beam 
polarized  in  the  x-direction,  is  obtained  by  setting  m  =  p  =  0  and  ej^  =  e^ 
in  Eq.  (14).  Using  Eq.  (14)  and  Eq.  (6),  the  Gaussian  beam  is  given  by 

E  =  E  exp(-r^/r^  +  i<|>)e  /2  +  c.c.,  (23) 

*“0  S  X 

where  =  k^  +  0  +  otr^/r^  and  the  functions  E  ,  9  ,  a,  and  r  are  given  by 

*  O  S  O  O  s 

Eqs.  (21)  and  (22)  with  p  =  0.  From  Eq.  (22a), 

EQ(^,q)  =  E^(^,q  =  0)rg(^,q  =  0)/rg(^,q),  where  E^(^,q  =  0)rg(^,q  =  0)  is 
proportional  to  the  optical  beam  power  P(^).  The  intensity  and  power 
associated  with  the  Gaussian  beam  in  a  medium  of  refractive  index  n^  are 
respectively 

I  =  (c/4n)<E  X  B>-e,  =  I  exp(-X),  (24a) 

~  ~  z  p 
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P  =  (n/2)r^  j  dX  KX)  =  (n/2)r^Ip,  (24b) 

o 

2 

where  Ip  =  (cn^/8ii)E^  is  the  peak  intensity  along  the  axis,  r  =  0  and  <  > 
denotes  time  averaging. 

To  determine  the  optical  beam  dynamics  from  Eqs.  (21)  and  (22)  the 
source  functions  G(h)  and  H(h)  in  Eq.  (20)  are  evaluated  with  p  =  0. 
Substituting  Eq.  (13)  into  Eq.  (20)  and  integrating  over  X  we  obtain 


=  -(k/2) 


l/2j 


-  (k  /k)' 
'  po 


"G1 


vOg2^-' 


(25) 


where  the  functions  Og2  represent  filling  factors  which  are 

essentially  the  ratio  of  the  cross-sectional  area  of  the  plasma  to  that  of 
the  optical  mode  times  the  normalized  plasma  density.  The  functions 
and  Og2  given  by 


‘^Gl 

=  f  dX  (k„/k^^)^ 

'  1  ' 

^"G2^ 

J  P  po 

0 

a-x> 

exp(-X). 


(26) 


Ionization  is  maximum  where  the  optical  field  amplitude  is  maximum,  i.e., 
at  r  =  0  for  the  fundamental  Gaussian  beam.  Since  the  tunneling  ionization 
rate  W(|e|)  depends  exponentially  on  the  field  amplitude,  the  radial 
profile  of  the  plasma  density  will  be  highly  peaked  about  the  axis  r  =  0. 
Equation  (26)  can  be  simplified  by  expanding  the  integrand  about  r  =  0, 
which  gives  =  0^2  =  where 

0  “ 

Og  =  J  d^'Kg(^')(r2(^')/r2(0]J  dX  exp  [-bg(  )X/2] 

^  0 


o 

'  d^'Kg(^')[r^(^')/rJ(0)/b(,(O, 


{27) 
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and 

Kg(0  =  4n^(3/n)^^^(S^/c)(Uj)^^'^[E^(0]'^^^exp[-bg(0],(28a) 

bg(0  =  (2/3)(Uj)^^^/E^(0,  (28b) 

and  E  =  E  /E„.  The  quantities  K^,  bp,  r^,  and  E  are  functions  of  ^  and 

O  O  H  b  b  S  o 

2  2 

ri,  whereas  the  power  P  ~  r  E  is  only  a  function  of  ^  as  implied  by  Eq. 

s  o 

(22a). 


A.  Envelope  Equation  for  Gaussian  Beam 

Using  Eqs.  (25)-(28),  the  envelope  equation  for  the  Gaussian  beam  in 
Eq.  (21)  becomes 


-  P/P^G  +  0/2)f^po‘'G)- 

2  2 

where  r^^  is  constant,  =  ^^^0/2  =  nn^rg^/X  is  the  Rayleigh  length 

associated  with  the  spot  size  r  ,  R  =  r  /r  is  the  normalized  spot  size, 

SO  s  so 

P  =  (n/2)I  r^  is  the  total  power  and  Pjjg  =  X^/( 2 101^02)  is  the  nonlinear 
focusing  power  for  the  Gaussian  beam.  The  terms  on  the  right-hand 

side  of  Eq.  (29)  denote,  respectively,  vacuum  diffraction,  nonlinear 
focusing  and  plasma  defocusing. 

In  the  absence  of  ionization  (o^  =  0)  the  envelope  equation  in  Eq. 
(29)  has  the  solution 


r  =  r 
s  so 


1/2 

(1  .  <1  -  P/PNoXt^ZEo)^]  • 


(30) 


where  r_  is  assumed  to  be  zero  at  h  =  0.  For  P  <  P^jq  the  optical  beam 


diffracts  with  an  effective  Rayleigh  length  given  by 


Zr  =  (1  -  P/PnG> 


(31) 
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For  P  =  diffractive  spreading  balances  nonlinear  focusing  and  a 

matched  self-guided  beam  can,  in  principle,  be  obtained.  However,  small 
changes  away  from  P  =  Pj^g  will  result  in  loss  of  equilibrium.  For  P  >  Pj^g 
the  optical  beam  self-focuses.  In  the  absence  of  ionization  the  beam 
focuses  down  to  zero  spot  size  with  a  focal  length  given  by 

However,  as  the  beam  focuses  the  intensity  on  axis  increases  resulting  in 
ionization  and  plasma  defocusing,  as  is  described  by  Eq.  (29). 


B.  Self-Guided  Gaussian  Beam 

In  the  presence  of  ionization,  self-guided  solutions  to  Eq.  (29)  can 

2  2 

be  obtained.  The  condition  for  a  self-guided  beam,  i.e.,  3  R/3lr|  =  0,  is 

P/P„_  -  1  =  klrla„/2  >  0.  (33) 

NG  po  S  G  ” 

Upon  taking  the  derivative  of  this  expression  with  respect  to  ^  and  using 
Eq.  (27)  we  find  that  for  a  self-guided  beam 

9P/85  .  -  lc2^P„gr2(0K5(0/bg(S) 

=  -  Phg(0,  (34) 


where 

>’g  - 

is  a  function  of  E^(0-  The  solution  of  Eq.  (34)  yields  the  self-guided 
optical  Gaussian  beam  power  as  a  function  of 


P(^)  =  PfjG^^P 


a  )dl 


I 


(35) 
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Equation  (33)  or  (35)  describes  a  family  of  equilibrium  solutions,  i.e., 

there  are  various  equilibrium  profiles  P(0>  which  satisfy 

these  equations.  For  example,  if  an  equilibrium  is  chosen  such  that  E^(^) 

is  constant  along  the  optical  beam,  then  h^  is  constant  and  Eq.  (35) 

2 

implies  P(^)  =  Pjjgexp(-yig)  and  the  spot  size  profile  is  given  by  = 

(16/cn^)P(0/E^-  Behind  the  beam  front,  ^  <  0,  the  optical  beam  power  and 
plasma  density  increase  such  that  the  nonlinear  self-focusing  term  and  the 
plasma  defocusing  term  remain  balanced.  Other  types  of  equilibria  can  be 
found,  for  example,  one  in  which  t^g(0  is  constant. 


C.  Stability  of  Gaussian  Beam 

In  this  section,  the  self-guided  beam  equilibrium  described  above  is 

shown  to  be  inherently  unstable,  i.e.,  the  beam  will  undergo  what  we  refer 

to  as  an  ionization-modulation  (IM)  instability.  The  IM  instability  is  due 

to  varying  degrees  of  ionization  along  the  beam  and  results  in  the 

modulation  of  the  beam  envelope  and  the  disruption  of  the  back  of  the  beam. 

To  examine  the  stability  of  the  self-guided  beam  equilibrium,  the  envelope 

equation,  Eq.  (29),  is  expanded  about  the  equilibrium  solution.  The 

perturbations  6r(^,  tl)  and  5E(^, h)  such  that  rg(^)  +  6r  and  E^(^)  +  SE 

denote  the  perturbed  spot  size  and  optical  field  amplitude  respectively. 

Furthermore,  since  the  optical  beam  power  within  the  paraxial  approximation 

is  nonevolving,  i.e.,  independent  of  Vi»  the  perturbations  Sr  and  SE  are 

related  by  Sr  =  -  (r  /E  )SE.  Expansion  of  the  envelope  equation,  Eq.  (29), 

s  o 

yields 


2  2 
a  Sr/ an  = 


4k^ 


-E2_ 


1  2  3,i.\ 
k  r^a) 


d^'Kg(5')rg(^')Sr(^',n), 


(36) 
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where  bg  »  1  has  been  assumed  (typically  the  case).  For  the  special  case 
of  an  equilibrium  with  a  constant  spot  size,  r  (^)  =  r  ,  Eq.  (36)  becomes 

S  50 


’  I  - 


L 


8r  =  0, 


(37) 


where  =  f2k  /(kr  )1  K_(0-  Equation  (37)  can  be  solved  taking  a 
g  V,  po  so  )  \j 

Laplace  transform  in  the  h  variable,  yielding 


Sr 


r  ®  1 

J  ds  exp 

sh  -  s"^J  d^  kg 

B 

K  ■* 

(38) 


where  s  is  the  Laplace  transform  variable  and  the  integration  is  over  the 
Bromwich  contour.  The  asymptotic  behavior  of  Sr  can  be  found  by 
integrating  Eq.  (38)  using  the  saddle  point  method. 


Sr  ~  exp 


I  (l  ±  1-J3](2r,2 

I 


1/3 


(39) 


Alternatively,  the  asymptotic  behavior  of  Sr  can  be  determined  from 

1/3  2/3 

Eq.  (37)  by  assuming  Sr  is  a  function  of  only  the  variable  x  =  (-0  h 
Substituting  Sr  =  Sr(x)  into  Eq.  (37)  yields 

(4/27)  [a^/ax^  +  (3/2)x”^a^/ax^  -(l/2)x"^a/axjsr  =  -kgSr-  (^0) 

In  the  asymptotic  limit,  x  Eq.  (40)  reduces  to 

[a^/ax^jSr  =  -(27/4)k^Sr,  (41) 

which  yields  the  solution  in  Eq.  (39). 

The  growth  rate  in  Eq.  (39)  can  be  simplified  by  noting  that,  along 
the  axis  r  =  0,  the  equilibrium  plasma  density  profile  is  given  by  an^/a^ 
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and  Eq.  (39)  becomes  6r  -  exp|^(l  ±  i>r3)N^(  ri)  j ,  where 

is  the  number  of  e-folds.  If  the  equilibrium  is  nearly  constant  in  the 
plasma  density  profile  is  given  by  n^/n^^  ^  I^I^G  number  of  e-folds 

is  where  j„ 

instability  grows  as  a  function  of  the  distance  behind  the  head  of  the 
optical  beam,  |^|,  and  the  propagation  distance  tj. 

The  dependence  of  on  ^  indicates  that  the  number  of  e-folds  at  the 

back  of  the  beam  is  greater  than  near  the  front.  The  IM  instability 

disrupts  the  back  of  the  beam,  and  the  disruption  point  propagates  toward 

the  front.  The  disruption  point  can  be  defined  as  the  point  on  the  beam 

where  the  initial  perturbation  is  increased  by  exp(N^),  where  »  1  is 

the  number  of  e-folds  necessary  for  disruption.  This  point  moves  toward 

the  front  of  the  beam  with  relative  velocity  =  “  c(  3N^/3tl)/(  3N^/3^) , 

where  =  N^.  For  the  case  where  the  plasma  density  profile  is 

linear,  i.e.,  n  =  K-n  Ul,  the  disruption  velocity  in  the  beam  frame  is 
p  u  no  '  ' 

V,  =  2cN^a”^n~^.  (^3) 

d  o  o 

To  gain  some  understanding  of  the  IM  instability,  consider  increasing 
the  spot  size  of  an  initially  matched  optical  beam,  i.e.,  5r(h  =  0)  >  0. 

In  this  case  the  beam  intensity  and  ionization  rate  are  reduced  resulting 
in  less  plasma  generation  and  enhanced  focusing  of  the  beam.  The  focusing 
optical  beam  overshoots  its  equilibrium  value  such  that  Sr  <  0  some 
distance  behind  the  beam  front.  When  Sr  <  0,  the  intensity,  ionization 
rate,  and  plasma  density  increase,  causing  the  beam  to  defocus  and 
overshoot  its  equilibrium  value.  This  focusing  and  defocusing  of  the  beam 
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results  in  the  IM  instability.  The  modulation  amplitude  and  period  are 
functions  of  the  distance  back  from  the  head  of  the  optical  beam,  |^|,  and 
the  propagation  distance,  Vl  =  z,  as  indicated  by  Eq.  (42). 

D.  Numerical  Results  for  Gaussian  Beams 
1.  Dynamic  Solutions 

The  propagation  dynamics  of  the  fundamental  Gaussian  beam  is  studied 
by  numerically  solving  the  envelope  equation,  Eq.  (29),  for  the  spot  size 
rg(^,D)«  The  envelope  equation  is  numerically  integrated  in  r\  (axial 
propagation  distance)  using  a  finite-difference  method,  where  the  initial 
conditions  rg(^,iri  =  0)  and  are  specified.  Note  that  in 

t 

evaluating  the  filling  factor  OQ(^,ri)  in  Eq.  (27),  the  integral  over  ^  is 
carried  out  at  every  ti  step,  since  the  integrand  in  Eq.  (27)  is  a  function 

r 

of  ^  and  h-  We  consider  a  linearly-polarized  laser  pulse  with  a  Gaussian 

radial  profile  and  an  initial  (D  =  z  =  0)  axial  profile  given  by  E^(^,0)  = 

E  ^sin(ii|^|/L)  for  -L  <  ^  <  0,  where  Ep^  =  (Slllp^/c)  is  the  initial  peak 

13  2 

electric  field,  I  =  3.0  x  10  W/cm  is  the  initial  peak  intensity,  and  L 
po 

=  60  ym  is  the  pulse  length.  With  wavelength  X  =  1  ym  and  initial  spot 

size  r^^  =  80  ym,  the  peak  power  is  =  3.0  GW  and  the  diffraction  length 

is  Z_  =2.0  cm.  The  optical  pulse  propagates  in  air  at  1  atm:  neutral 
KO 

19  -3  37  -19 

gas  density  n^^  =  2.7  x  10  cm  ,  nonlinear  index  n2  =  5.6  x  10 

cm^/W,  normalized  ionization  potential  Uj/Ug  =  1.07,  and  nonlinear  focusing 

power  Pjjg  =  2.8  GW  (Pq/P^G  “ 

The  simulation  begins  with  the  optical  pulse  at  focus  (Sr^/Sti  =  0)  in 
the  neutral  gas.  With  the  initial  value  of  the  filling  factor  computed 
via  Eq.  (27),  the  envelope  equation,  Eq.  (29),  is  integrated  in  the 
simulation  variables  ^  =  z  -  vt  and  q  =  z.  Figures  2(a)  and  2(b)  show  the 
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initial  h  =  0  optical  beam  intensity  I  and  plasma  density  n^  versus  (r,^). 

In  Figs.  2(a)  and  2(b),  the  direction  of  propagation  is  towards  the  right. 

Plots  of  I  and  n  versus  radius  at  the  pulse  center  (^  =  -  30  pm)  are  shown 
P 

in  Fig.  2(c)  for  this  case.  Note  also  that  the  nonlinear  nature  of  the 
ionization  process  causes  the  plasma  density  gradient  versus  both  r  and  I, 
to  be  considerably  sharper  than  the  intensity  gradient. 

The  evolution  of  the  optical  pulse  is  shown  in  Fig.  3(a-d),  where  the 
spot  size  r^  (dashed  line),  intensity  I  on  axis  (solid  line),  and  plasma 
density  np  on  axis  (dotted  line)  are  plotted  versus  ^  at  (a)  h  =  2  =  0,  (b) 
6  cm,  (c)  8  cm,  and  (d)  10  cm.  Initially,  the  spot  size  is  constant  along 
the  optical  pulse,  as  shown  in  Fig.  3(a).  Because  >  Pj^g,  the  center  of 
the  pulse  is  focused  while  the  front  and  back  portions  diffract,  as  seen  in 
Fig.  3(b).  At  ^  =  -25  pm,  where  P  *  Pj^g  and  n^  «  n^^,  diffraction 
balances  nonlinear  focusing  and  the  spot  size  remains  constant  at  r^  =  r^^. 
Behind  this  point,  focusing  increases  the  optical  intensity,  producing  a 
corresponding  increase  in  the  ionization  rate.  Because  ionization  is  a 
highly  nonlinear  process,  the  steepness  of  the  plasma  density  gradient  also 
increases.  Increased  ionization  and  increased  plasma  density  gradients  are 
shown  in  Figs.  3(b-d).  Increased  ionization  causes  the  latter  portion  of 
the  optical  pulse  to  diffract,  as  can  be  seen  in  Figs.  3(c,d).  The  rapid 
change  in  the  plasma  density  at  the  steepening  ionization  front  results  in 
a  correspondingly  rapid  change  in  the  focusing  of  the  optical  pulse.  This 
results  in  an  increasingly  narrow  intensity  spike  at  the  ionization  front. 

The  optical  pulse  structure  observed  in  Figs.  3(c,d)  occurs  even  when 
the  power  P  greatly  exceeds  the  nonlinear  focusing  threshold  Pj^g.  For 
example.  Figs.  4(a-d)  show  the  evolution  of  a  pulse  with  =  2. 

Except  for  the  initial  power  and  the  initial  spot  size  r^^  =  110  ym  (and  a 
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corresponding  change  in  =  3.8  cm),  the  parameters  of  Fig.  4  are 

identical  to  those  of  Fig.  3.  As  in  Fig.  3,  there  is  a  point  near  the 

front  of  the  pulse  in  Fig.  4  (at  ^  =  -15  pm)  where  P  =  n^  «  n^^,  and 

diffraction  balances  nonlinear  focusing.  Also  as  before,  nonlinear 
focusing  of  the  pulse  behind  this  point  leads  to  an  increasingly  steep 
ionization  front.  This,  in  turn,  produces  an  increasingly  narrow  optical 
intensity  spike. 

2.  Self-Guided  Solutions 

Examples  of  matched  beam  equilibria  are  shown  in  Figs.  5  and  6.  In 

both  cases,  we  consider  a  linearly  polarized  X  =  1  pm  optical  pulse  with  a 

19  -3 

Gaussian  radial  profile  propagating  in  air  at  1  atm  (n^^  =  2.7  x  10  cm  , 

=  5.6  xlO^^  cm^/W,  Uj/Ujj  =  1.07,  and  Pj^g  =  2.8  GW).  Figure  5  shows 

optical  power  profiles  (solid  lines)  and  plasma  density  profiles  (dashed 

lines)  plotted  versus  ^  along  the  axis  for  equilibria  with  constant 

profiles.  Equilibria  are  shown  for  three  different  values  of  the  optical 

intensity:  1^  =  =  5  x  10^^  W/cm^,  I2  =  6  x  10^^  W/cm^,  and  =  7  x 

10^^  W/cm^.  Here,  =  (8JlIp/c)^'^^  is  constant  along  the  length  of  the 

optical  pulse,  such  that  the  variation  in  power  P(^)  corresponds  to  a 

1/2 

variation  in  spot  size  r  =  (2P/liI  )  .  Note  that  the  constant  E  profile 

s  p  ^ 

produces  a  constant  ionization  rate  and  a  linear  rise  in  Op/n^^  «  1. 

Also,  the  power  profiles  are  exponential  functions  as  given  in  Eq.  (35). 

Figure  6  shows  optical  power  and  plasma  density  profiles  for 

equilibria  with  constant  =  r^^  profiles.  In  this  case,  matched  E^(0 

profiles  are  determined  numerically  from  Eq.  (34)  for  three  different 

13  2 

values  of  the  leading-edge  (^  =  0)  intensity:  Ip  =  =  5  x  10  W/cm  ,  I2 

=  5.1  X  10^^  W/cm^,  and  I^  =  5.2  x  10^^  W/cm^.  In  this  case,  the  variation 
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in  optical  power  P(J;)  corresponds  to  a  variation  in  intensity  Ip  =  2P/iir^^, 

such  that  E  increases  with  k|  along  the  length  of  the  pulse.  As  a 
o 

result,  the  ionization  rate  increases  as  a  function  of  |^|.  Increased 

ionization  (defocusing)  requires  increased  power  (focusing)  to  compensate, 

further  increasing  the  ionization  in  a  highly  nonlinear  manner.  As  a 

result,  the  constant-r  equilibrium  profiles  can  be  very  sensitive  to  the 

s 

value  of  Ip  as  in  Fig.  6. 

3.  Ionization-Modulation  Instability 

An  example  of  the  IM  instability  for  a  fundamental  Gaussian  beam 

obtained  by  numerical  solution  of  the  envelope  equation,  Eq.  (29),  is  shown 

in  Fig.  7.  Here,  we  consider  the  propagation  in  air  of  a  constant-r^ 

equilibrium  with  Ip(^  »=  0)  =  3.0  x  10^^  W/cm^,  r^  =  r^^  =  78  ym,  and  = 

1.9  cm.  In  this  case,  there  is  very  little  initial  ionization  and  the 

growth  of  the  instability  is  extremely  slow  with  P(0  “  Pjjq  =  2.8  GW  along 

the  length  of  the  optical  pulse.  The  evolution  of  the  optical  pulse  is 

shown  in  Figs.  7(a-f),  where  the  spot  size  r^  (solid  line)  and  plasma 

density  n  on  axis  (dashed  line)  are  plotted  versus  ^  at  (a)  n  =  z  =  0,  (b) 
P 

400  cm,  (c)  450  cm,  (d)  500  cm,  (e)  550  cm,  and  (f)  600  cm.  In  Fig.  7,  the 
direction  of  propagation  is  towards  the  right. 

The  simulation  begins.  Fig.  7(a),  with  the  optical  pulse  at  focus 
(9r^/3Vl  =  0)  in  the  neutral  gas.  In  Fig.  7(a),  the  spot  size  r^  is 
constant  along  the  pulse  and  np(0  increases  linearly  since  E^  is 
approximately  constant.  At  later  times.  Figs.  7(b-d),  oscillations  in  r^ 
cause  oscillations  in  the  ionization  rate  such  that  each  region  where  r^ 
has  decreased  corresponds  to  an  increase  in  ionization.  This  is 
particularly  noticeable  at  the  back  of  the  pulse  (^  =  -60  ym)  in  Fig.  7(c). 
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Eventually,  there  is  a  large  enough  increase  in  the  plasma  density  so 

that  the  latter  portion  of  the  optical  pulse  is  defocused,  i.e.,  the 

guiding  is  disrupted.  When  the  optical  pulse  is  sufficiently  defocused  the 

ionization  rate  falls  and  dn^/dJI  =  0.  Thus,  an  "ionization  front"  develops 

which  propagates  forward  in  the  beam  frame.  This  can  be  seen  in  Fig.  7(e), 

where  the  ionization  front  is  at  t  =  -40  pm,  and  in  Fig.  7(f),  where  the 

ionization  front  is  at  ^  -  -33  pm.  Figures  7  indicate  that  the  disruption 

velocity  is  in  good  agreement  with  Eq.  (43). 

The  growth  of  the  instability  of  Fig.  7  is  plotted  versus  I  at  fixed  q 

=  z  =  550  cm  in  Fig.  8(a),  where  ln|Ar|),  from  the  numerical  integration  of 

the  envelope  equation,  is  compared  to  the  number  of  e-folds  from  Eq. 

(42).  Here,  Ar  =  (r  -  r  )/r  .  Similarly,  ln|ar|)  versus  q  =  z  at  fixed 

s  so  so 

^  =  -49  pm  is  shown  in  Fig.  8(b).  For  both  plots,  excellent  agreement  is 
observed  between  the  slope  of  and  the  peaks  of  the  ln|Ar|)  curve.  As 

expected,  agreement  tends  to  break  down  for  q  -^0,  where  the  growth  is 
not  yet  asymptotic,  and  for  ln|Ar|)  -*  0,  where  the  growth  is  nonlinear. 

IV.  Higher-Order  Radially-Polarized  Beam  Propagation 

We  now  consider  the  dynamics  of  a  radially-polarized  optical  beam 
propagating  in  a  gas  undergoing  ionization.  The  radially-polarized  optical 
beam  is  formed  by  taking  m  =  0,  p  =  1  in  Eq.  (14),  setting  aQ  =  bp  = 

E  exp(i6  ),  e.  =  e  and  e!  =  e  .  The  resulting  field,  from  Eqs.  (6)  and 

0  0  -L  X  JL  y 

(14),  is 

E  =  E^(-f2r/rg)exp(-  +  i'l')ej^/2  +  c.c.,  (44) 
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where  +  6  +  ar^/r^  and  e^  =  cosG  4  +  sin©  e  is  the  unit  radial 

o  s  r  X  y 

vector.  The  functions  E  ,  0  ,  a,  and  r  satisfy  Eqs.  (21)  and  (22)  with  p 

0  0  s 

=  1.  From  Eq.  (22a),  E^(^,n)  =  =  0)/rg(£;,n),  where 

E^(^,ri  =  0)r^(^, n  =  0)  is  proportional  to  the  optical  beam  power  P(0-  The 
o  s 

axial  field  component  associated  with  the  radial  field  in  Eq.  (44),  as 
obtained  from  V*E  =  0,  is  maximum  along  the  axis  and  given  by 

E^  -  Eq[i  -  (r^/r^)(l  -  ia)jexp(-  +  i^)e^/2  +  c.c., 

(45a) 

E  (r  =  0)  *  ^  E^exp(ik^  +  i0^)e^/2  +  c.c.,  (45b) 

s 

where  the  expression  in  Eq.  (45b)  is  valid  along  the  z-axis,  i.e.,  r  =  0. 
The  intensity  and  power  associated  with  the  radially-polarized  beam  in  a 
medium  of  refractive  index  n^  are  respectively 

I  =  (c/4ll)<E  X  B>-e,  =  I  Xexp(l-X),  (46a) 

^  z  p 

00 

P  =  (n/2)rg  J  dx  I(X)  =  (eii/2)rglp,  (46b) 

0  . 
n 

where  I  =  (cn  /8n)E'^/e  is  the  maximum  intensity,  which  occurs  at  X  =  1 
P  o  o 

(r  =  r  /-JT),  e  =  2.72  and  <  >  denotes  time  averaging.  Similar  higher-order 
s 

radially-polarized  modes  have  been  produced  using  an  axicon  focusing 

..  38,39 

configuration. 

To  determine  the  optical  beam  dynamics  from  Eqs.  (21)  and  (22)  the 
source  functions  G(ri)  and  H(Yi)  in  Eq.  (20)  are  evaluated  with  p  =  1. 
Substituting  from  Eq.  (13)  into  Eq.  (20)  and  integrating  over  q  yields 

/r*  \  r 

G  J-  9  ^ 

.  -(k/2)  (e/axnj/n^)!  .  ,  (47) 
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where  the  filling  factors  2  given  by 


"R1 


^<'r2 


J  « 

0 


'  1  ' 

a-x/2> 


(48) 


The  tunneling  ionization  rate  W(  |e|)  is  maximum  when  the  optical  field  is 
maximum.  The  magnitude  of  the  radial  polarized  field  in  the  expression  for 
kp,  Eq.  (12),  can  be  written  as  |e|  =  E^(l  +  y)  exp[(l  +  y)/2],  where  y  = 
2r^/r^(^  )  -  If  and  has  a  maximum  at  y  =  0,  i.e.,  r  =  r  (^  )/’f2.  Expanding 

S  s 

III  for  y  «  1  gives  |e|  =  (E^/'fe)(l  -  y  /4).  Using  this  expression  in 
W(|e|),  and  noting  that  W(|e|)  is  highly  peaked  about  y  =  0,  gives 


2  2  r°  '  '  -b  (5')y2/4 

kp(0/kpo  =  ej  d^  Kj^(^  )e  ^  ,  (49) 


where 


.-1/2 


K^(0  =  4n^(3/n)^^V^^‘^(Q^/c)(0j)^^‘^(E^(O]  exp[-bj^(0), 

(50a) 


bR(0  =  I  e^^^(Uj)^^^/E^(0,  (50b) 

and  E  =  E  /E„.  Inserting  Eq.  (49)  into  the  expressions  for  the  filling 

O  On 

factors,  Eq.  (48),  and  assuming  b^  »  1  allows  the  integration  over  r  to  be 
carried  out  yielding 


R1 


'.0R2/ 


=  2fi  J  d^'KR(^')b‘^''^(^') 


U-X^/2>l 

s 


Xgexp(l  -  Xg),  (51) 


where  X^  =  )/r^(K)- 
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A.  Envelope  Equation  for  Higher-Order  Beam 

Using  Eq.  (47)  for  the  source  function  H,  the  envelope  equation  in  Eq. 
(21)  becomes 

ih/ix?-  .  Z-2r-3(i  -  P/P^  .  (52) 

2  2 

where  r  is  constant,  Z-  =  kr  /2  =  nn  r  /X  is  the  Rayleigh  length 
so  Ro  so  o  so 

associated  with  the  spot  size  r  ,  R  =  r„/r^^  is  the  normalized  spot  size, 

^  SO  s  so 

2  2 

P  =  (en/2)Iprg  is  the  total  power,  and  Pj^  =  APj^g  =  2X  /(nn^n2)  is  the 

nonlinear  focusing  power  for  the  higher-order  radially-polarized  beam.  The 
terms  on  the  right-hand  side  of  Eq.  (52)  denote  respectively,  vacuum 
diffraction,  nonlinear  focusing  and  plasma  defocusing.  In  the  absence  of 
ionization,  aj^2  =  solution  to  Eq.  (52)  for  r^  is  given  by  Eq.  (30) 

with  Pjjg  replaced  by  Pj^. 

B.  Self-Guided  Higher-Order  Beam 

2  2 

The  condition  for  a  matched  beam,  i.e.,  3  R/3V1  =  0,  is  given  by 

f(fNR  -  1  >  0-  (53) 

Equation  (53)  describes  a  family  of  equilibrium  solutions,  i.e.,  there  are 

various  equilibrium  profiles  P(^),  E  (0>  and  r  (^)  which  satisfy  Eq.  (53). 

o  s 

For  the  special  case  of  a  constant  spot  size  matched  beam,  r^  =  the 

matched-beam  power  is  given  by 

3P/8C  .  -(-n^2)k2^P^4l4,(5)/bl''2(j;,,  (54) 

where  1^(0  and  are  functions  of  P(0  through  E^(0»  i-e.,  P(0  = 

(cno/16)E2(t)4. 
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C.  Stability  of  Higher-Order  Beam 

The  stability  of  the  matched-beam  equilibrium  can  be  analyzed  by 
perturbing  the  envelope  equation^  Eq.  (52),  about  the  matched  beam 
solution.  This  is  accomplished  by  introducing  a  constant  power 
perturbation  to  the  spot  size  and  field  amplitude  of  the  form  r^(^)  + 
8r(^,ri)  and  +  SE(^,ri),  where  SE/E^  =  -Sr/r^.  In  the  limit  bg  »  1, 

the  envelope  perturbation  satisfies  the  equation 


2  2 

3  8r/3ri  = 


1  2 

f—  K  r 

sn  po  so 
2  2  3 

Ro 


(2  - 


Xpbg(^  )Sr(0 


•  2X3(2 


4X2  ^ 


1-X 

-1/2.. '.,,2/  *s 


>'>r  >2;* 


(55) 


For  a  constant  spot  size  equilibrium,  X  =  1,  Eq.  (55)  becomes 

s 


35r 

3^ 


=  kJCOSr, 


(56) 


where  =  2'rn[kp^/(krg^)]  In  the  limit  bg  »  »  1,  where 

N  is  the  number  of  e-folds  of  the  instability,  the  second  term  on  the  left 
e 

of  Eq.  (56)  can  be  neglected  and  the  asymptotic  behavior  of  Sr  is  given  by 


Sr  ~  exp 


I  (1  1  i-ls)  (2(1 


2 


(57) 


Notice  that,  for  an  equilibrium  in  which  E^(0  is  nearly  constant,  the  peak 
equilibrium  plasma  density  occurs  at  r  =  r^/i2  and  is  given  by  - 
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-eKj^^.  In  this  limit,  the  asymptotic  behavior  of  6r  can  be  written  as  Sr  ~ 
exp^(l  ±  i-fT)Ngj,  where 


,1/3 


(n  /n  )  (k  r  h/Zp 
'  p  no'  '  po  so  Ro 


(58) 


is  the  number  of  e-folds.  Hence,  the  growth  rate  of  the  IM  instability  for 

2  1/6 

the  higher-order  optical  beam  differs  roughly  by  a  factor  of  (iibj^/4e  ) 

1  from  that  of  the  fundamental,  assuming  equal  values  at  Op/n^^^, 

and 


D.  Numerical  Results  for  Higher-Order  Beam 
1.  Dynamic  Solutions 

The  propagation  dynamics  of  a  higher-order  radially-polarized  laser 

pulse  can  be  described  by  a  numerical  solution  of  the  envelope  equation, 

Eq.  (52).  Initially,  n  =  z  =  0,  the  profile  of  the  radial  electric  field 

is  |E(r,0|  =  E^(^,0)(42r/rg)exp(-r^/r^),  with  an  axial  profile  E^(0  = 

E  e^^^sin(nk|/L)  for  -L  <  ^  <  0,  where  L  =  60  pm  is  the  optical  pulse 
po  '  '  —  — 

length.  The  peak  initial  field,  occurring  at  r  =  is  Ep^  = 

(Sul  /c)^'^^,  where  I  =  4.7  x  10^^  W/cm^  is  the  peak  initial  optical 
po  po 

intensity.  With  wavelength  X  =  1  pm  and  initial  spot  size  r^^  =  35  pm,,  the 
peak  power  is  =  2.5  GW  and  the  diffraction  length  is  =  0.4  cm.  We 
consider  propagation  in  Hydrogen  (H2)  at  30  atm:  nonlinear  index  n2  =  3.3 
X  10"^®  cm^/W,  ionization  energy  Uj  =  15.4  eV,  neutral  density  n^^  =  8.1  x 
10^^  cm~^,  and  nonlinear  focusing  power  =  1.9  GW  (Pq/Pjjj^  “  1.3). 
Hydrogen  is  chosen  for  its  low  atomic  number  and,  hence,  Bremsstrahlung 
losses  can  be  neglected,  see  Appendix  C. 
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The  simulation  begins  with  the  optical  pulse  at  focus  (9r  /Sh  =  0)  in 

s 

the  neutral  gas  with  intensity  I(r,?;)  and  plasma  density  np(r,0  profiles 

initialized  as  shown  in  Fig.  9(a,b),  where  the  direction  of  propagation  is 

towards  the  right.  The  accelerating  field  =  |E^|  profile  is  plotted  in 

Fig.  9(c),  where  E  =  E  exp(ikO/2  +  c.c.  with  E  given  by  Eq.  (45).  Plots 

z  z  z 

of  I,  E^,  and  np  versus  radius  at  the  pulse  center  (^  =  -30  ym)  are  shown 
in  Fig.  9(d)  for  this  case.  The  highly  nonlinear  nature  of  the  ionization 
process  causes  the  plasma  density  gradient  versus  both  r  and  K  to  be 
considerably  sharper  than  the  intensity  gradient. 

The  evolution  of  the  optical  pulse  is  shown  in  Figs.  lO(a-d),  where 

the  spot  size  r  (dashed  line),  accelerating  field  E  on  axis  (solid  line), 

s  z 

and  plasma  density  ttp  (dotted  line)  are  plotted  versus  ^  at  (a)  q  =  z  =  0, 

(b)  0.4  cm,  (c)  0.7  cm,  and  (d)  1.0  cm.  In  Figs.  lO(a-d),  np(^)  is 

evaluated  at  r  =  r  .  /-il,  where  r  .  is  the  minimum  value  of  r  (O. 

min  min  s 

Initially,  the  spot  size  is  constant  along  the  optical  pulse,  and  the  peak 

accelerating  field  is  E^  =  400  MV/m,  as  in  Fig.  10(a).  Since  the 

center  of  the  pulse  is  focused  while  the  front  and  back  portions  diffract, 

as  shown  in  Fig.  10(b).  The  increasing  optical  intensity  in  the  center  of 

-  1/2 

the  pulse  increases  the  accelerating  field  since  E^  ~  I  .  Diffraction 

balances  nonlinear  focusing  at  ^  =  -21  ym,  where  P  «  Pj^,  and  the  spot  size 

remains  constant,  r  =  r  .  Behind  this  point,  focusing  increases  the 

s  so 

optical  intensity,  producing  a  corresponding  increase  in  the  ionization 
rate.  Increased  ionization  and  increased  plasma  density  gradients  are 
shown  in  Figs.  lO(a-c).  Increased  ionization  causes  the  optical  pulse  to 
diffract,  reducing  E^  in  the  latter  portion  of  the  pulse,  as  shown  in  Figs. 
10(c,d).  The  rapid  change  in  the  plasma  density  at  the  steepening 
ionization  front  results  in  a  corresponding  rapid  change  in  the  focusing  of 
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the  pulse  which  results  in  an  increasingly  narrow  spike  in  at  the 

ionization  front.  The  peak  accelerating  field  is  E  =  1.1  GV/m  at  h  =  z  = 

z 

I. 0  cm  e  2.5Zj^^. 

2.  Self-Guided  Solutions 

An  example  of  a  self-guided,  higher-order  radially-polarized  beam  of 
wavelength  X  =  1  ym  propagating  in  Hydrogen  (H2)  at  30  atm  is  shown  in  Fig. 

II.  For  H2  at  30  atm,  n^^  =  8.1  x  10^^  cm  n2  =  3.3  x  10  cm^/W, 

=  1.1,  and  P.,T,  ®  1*9  GW.  The  matched  beam  conditions  can  be  determined 

from  Eq.  (53)  together  with  Eq.  (51)  for  a  given  axial  intensity  profile. 

For  equilibria  with  constant  E^(^)  profiles,  i.e.,  constant  peak  intensity 
10  2 

of  I  =  4.7  X  10  W/cm  ,  the  matched  profiles  for  power  P  and  plasma 
P 

density  n  at  r  =  versus  ^  are  shown  in  Fig.  11.  For  these 

-4 

parameters,  the  degree  of  ionization  is  small,  ‘  spot 

size  is  given  by  rg(^)  =  2P( ^)/(neIp) ,  which  gives  rg(^  =  0)  =  r^^^  =  31  ym 

at  the  front  of  the  beam.  The  on-axis  accelerating  field,  also  shown  in 

Fig.  11,  has  the  maximum  value  450  MV/m  at  the  front  of  the  optical  beam. 

Since  E  ~  nt  ~  n  ,  the  accelerating  gradient  can  be  increased  by 
z  z  no 

increasing  the  gas  pressure. 

3.  Ionization-Modulation  Instability 

As  discussed  in  Sec.  IV. C.,  the  optical  beam  undergoes  an  IM 
instability.  Numerical  simulations  of  Eq.  (52)  show  that,  with  a  1% 
initial  perturbation  of  the  spot  size,  the  IM  instability  significantly 
disrupts  the  equilibrium  beam  profile  of  Fig.  11  after  -  10  cm  of 
propagation.  The  growth  rate  of  the  IM  instability  is  a  highly  nonlinear 
function  of  the  optical  intensity  through  the  plasma  density.  Reducing  the 
intensity  in  the  example  of  Fig.  11  to  3.2  x  10^^  W/cm^  (which  reduces  E^ 
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to  300  MeV/tn)  results  in  a  matched  beam  with  very  little  ionization  (the 

plasma  density  is  reduced  by  a  factor  of  ~  130)  and  little  variation  in 

power  and  spot  size  along  the  length  of  the  optical  pulse.  Simulations 

show  that  with  a  IX  perturbation  in  the  spot  size,  the  pulse  propagates  >  1 

meter  without  significant  disruption.  The  behavior  of  the  IM  instability 

is  shown  in  the  following  numerical  example.  In  this  example,  however,  the 

instability  is  allowed  to  grow  from  numerical  noise  to  facilitate 

comparison  to  the  theoretical  growth  rate. 

Propagation  of  a  guided  pulse  in  H2  at  30  atm  is  shown  in  Fig.  12. 

13  2 

Initially,  Ip  =  3.2  x  10  W/cm  is  constant  throughout  the  pulse,  the  peak 
accelerating  field  is  =  300  MV/m,  the  spot  size  is  r^  «  37  ym  (Zj^^  * 

0.44  cm),  and  P(0  =  =  1*9  GV.  The  evolution  of  the  optical  pulse  is 

shown  in  Figs.  12(a-d),  where  the  spot  size  r  (dashed  line),  accelerating 

s 

field  on  axis  (solid  line),  and  plasma  density  np  at  r  =  (dotted 

line)  are  plotted  versus  t  at  (a)  h  =  z  =  0,  (b)  60  cm,  (c)  75  cm,  and  (d) 
90  cm.  Initially,  the  spot  size  is  constant  along  the  optical  pulse,  and 
the  peak  accelerating  field  is  E  =  300  MV/m,  as  in  Fig.  12(a).  At  later 

times.  Figs.  12(b-d),  oscillations  in  r  cause  oscillations  in  the 

s 

ionization  rate  such  that  each  region  where  r  has  decreased  corresponds  to 

s 

an  increase  in  -dOp/d^.  Eventually,  the  plasma  density  is  sufficiently 
large  so  that  the  latter  portion  of  the  optical  pulse  is  defocused.  In  the 
region  where  the  pulse  is  sufficiently  defocused  the  ionization  rate  falls, 
dnp/d^  =  0,  and  an  ionization  front  develops  which  propagates  toward  the 
front  of  the  beam  at  the  disruption  velocity  v^,  see  Eq.  (43).  This  can  be 
seen  in  Fig.  12(c),  where  the  ionization  front  is  at  ^  =  -42  ym,  and  in 
Fig.  12(d),  where  the  ionization  front  is  at  ^  ss  -29  ym. 
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The  growth  of  the  perturbed  radius  due  to  the  instability  is  plotted 
versus  ^  at  fixed  h  =  75  cm  in  Fig.  13(a),  where  ln|Ar|),  from  the 
integration  of  the  envelope  equation,  Eq.  (52),  is  compared  to  the  number 
of  e-folds  N  from  Eq.  (58).  Here,  Z^r  =  (r  -  r  )/r  ,  where  r  (^)  is 

Q  S  Sv  SU  o 

plotted  in  Fig.  12(c).  Similarly,  growth  versus  q  =  z  at  fixed  I,  =  -40  ym 
is  shown  in  Fig.  13(b).  For  both  plots,  good  agreement  is  obtained  between 
the  slope  of  N^(0  and  the  peaks  of  the  ln|Ar|)  curve.  As  expected, 
agreement  tends  to  break  down  for  q  -»  0,  where  the  growth  is  not  yet  in 
the  asymptotic  limit,  and  for  ln|Zir|)  -♦  0,  where  the  growth  is  nonlinear. 

V.  Self-Guided  Inverse  Cherenkov  Accelerator 

A.  Electron  Energy  Gain  in  Vacuum 

The  axial  field  component  of  a  higher-order  radially-polarized  optical 

beam,  which  has  an  amplitude  on  axis  of  E^^  =  (2-r2/krg)E^,  can  be  used  to 

9  40-43 

accelerate  an  injected  electron  beam  propagating  along  the  z-axis.  ’ 

In  vacuum,  the  phase  velocity  Vp^^  of  the  optical  beam  is  greater  than  c  and 
near  the  focal  point  is  Vpj^/c  a  1  +  2/kZj^^.  Since  Vpj^  >  c,  phase  slippage 
between  the  electrons  and  optical  beam  will  occur.  For  a  highly 
relativistic  injected  electron  in  vacuum,  the  slippage  distance  is 
defined  as  the  distance  over  which  the  electron  phase  slips  by  one-half  an 
optical  period,  Lg(Vpj^  -  c)/c  a  X/2,  which  gives  a  iiZj^^/2.  It  can  be 
shown  that  a  highly  relativistic  electron  interacting  with  the  axial 
optical  field  E  in  vacuum,  Eq.  (59),  from  z  =  -“  to  z  =  »  experiences  zero 

net  energy  gain.  This  result  is  a  particular  case  of  the  Lawson-Woodward 
44,45 

theorem. 

It  can  also  be  shown, however,  that  if  the  interaction  distance 
is  limited  by  placing  a  mirror  approximately  one  Rayleigh  length  from  the 
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focal  point,  i.e.,  -Zj^^  1  z  <  ®,  then  the  maximum  energy  gain  for  a  highly 

1/2 

relativistic  electron  is  W^[MeV]  =  P  [GW].  Material  damage 

considerations,  associated  with  the  introduction  of  an  optical  component 

9  42 

near  focus,  place  serious  limits  on  the  optical  intensity.  ' 

B.  Electron  Energy  Gain  in  Conventional  ICA 

By  introducing  a  neutral  gas  into  the  interaction  region,  as  in  the 

inverse  Cherenkov  accelerator  (ICA),  the  phase  velocity  of  the  optical  beam 

6—8 

can  be  reduced  and  phase  slippage  reduced.  In  the  conventional  ICA,  ~  the 

optical  beam  diffracts  and  the  effects  of  nonlinear  self-focusing  and 

ionization  are  neglected.  For  a  higher-order  Gaussian  beam  the  phase 

velocity  near  focus  is  Vpj^/c  «  1  +  2/kZj^^  -  (n^  -  l)/n^.  Typically  n^  -  1 

«  1  and  is  proportional  to  the  neutral  gas  density  n^^^.  Proper  choice  of 

n  can  result  in  v  .  <  c  and  the  reduction  of  phase  slippage,  however, 
no  pn  — 

diffraction  remains  an  important  limitation. 

The  energy  gain  in  a  conventional  ICA  driven  by  a  higher-order 
Gaussian  beam  can  be  calculated  as  follows.  Assuming  vacuum  diffraction 
and  neglecting  nonlinear  and  plasma  effects,  the  axial  electric  field  along 
r  =  0  in  the  conventional  ICA  is  given  by 

B,(r  .  0)  .  -  E^^(l  .  (5,^ 

where  E  =  (2-f2/kr  )E  (z  =  0)  is  the  peak  axial  electric  field, 
zo  so  o 

«|/  =  kz  -  ckt/n  -  2tan“^z/Z„  +  tj;  is  the  phase  and  v|;^  is  a  constant. 

^  o  Ro  o  o 

Equation  (59)  follows  from  Eq.  (45)  with  E^(J;,ri)  =  E^(^,tl  =  0)rg^/rg(r|) , 

2  2  1/2 

where  r^  =  ^^^(1  +  h  /Zj^)  and  t|  =  Z.  The  energy  gain  of  a  highly 
relativistic  electron  interacting  with  the  optical  field,  given  by  Eq. 

(59),  in  a  medium  with  linear  refractive  index  n^  is  given  by 
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00 


(60) 


Wg  =  q  J  dz  E^(r  =  0)  =  -  2nqEj,^AkZ^^exp(-AkZj^^) , 

—  00 

for  Ak  >  0  and  W  =0  for  Ak  <  0,  where  Ak  =  k(n^  -  l/p  )/n  =  k(n  -  1  - 

l/2y^)/n^,  Yq  =  (1  -  is  the  electron  relativistic  mass  factor,  = 

V  /c,  V  is  the  axial  electron  velocity  and  an  initial  phase  of  sin<</  =  1 

was  assumed.  In  obtaining  Eq.  (60)  the  electron  trajectory  was  taken  along 

2 

r  =  0  from  z  =  -  “  to  z  =  »,  where  z  =  g^ct  a  ct(l  -  l/2r^).  Equation  (60) 
is  valid  as  long  as  the  energy  gain  is  less  than  the  initial  energy,  « 
(y^  -  l)mc  .  For  highly  relativistic  electrons  in  vacuum  (g^  =  1  and  n^  = 
1),  Ak  vanishes  and  the  energy  change  is  zero,  =  0,  in  accordance  with 
the  Lawson-Woodward  theorem.  In  a  gas,  the  energy  gain  is  maximum  when 
AkZj^^  =  1,  i.e.,  when  n^  -  1  =  n^/kZj^^  +  1/2Yq  =  (1/2)  ^(X/iirg)^/n^  +  l/r^j- 
In  this  case,  the  maximum  energy  gain  in  the  conventional  ICA  is  given  by 

"max  -  -  <»> 

The  maximum  energy  gain  is  the  product  of  the  peak  axial  field  E^^  and 

(2n/e)Z,j  ,  which  is  the  effective  acceleration  length.  In  terms  of  the 
KO 

optical  power,  Eq.  (61)  can  be  written  as 

W  [MeV]  =  2.3  P^^^[GW].  (62) 

max‘  ^ 

A  similar  result  can  be  obtained  in  vacuum  acceleration  by  limiting  the 
interaction  region  to  approximately  2Zj^  about  the  laser  focus  through  the 
use  of  optical  components  (e.g.,  mirrors). 

In  addition  to  the  higher-order  Gaussian  optical  beam  discussed  above, 
a  conventional  ICA  could  be  driven  by  a  first  order  Bessel  (axicon) 
beam.^’®  A  nonideal  Bessel  beam  (finite  transverse  extent)  consisting  of  N 
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2 

rings  (lobes)  will  propagate  a  distance  =  ^^Ro’  ^Ro  ” 

r^  is  the  width  of  an  individual  ring.  Consequently,  the  maximum  energy 

gain  in  a  conventional  ICA  driven  by  a  Bessel  beam  is  approximately  a 
1/2 

factor  of  N  greater  than  an  ICA  driven  by  a  higher-order  Gaussian  beam, 
assuming  equal  total  power  in  each  of  the  beams,  as  discussed  in  Appendix 
D. 

The  above  calculation  of  the  energy  gain  in  the  conventional  ICA 
assumes  vacuum  diffraction,  i.e.,  the  effects  of  nonlinear  self-focusing 
and  ionization  are  neglected.  Equation  (60)  only  holds  for  powers  below 
the  nonlinear  focusing  power,  P  «  Pj^,  and  intensities  below  the 
ionization  threshold,  I  «  Ij.  Typically,  at  atmospheric  pressures,  Pj^j^  ~ 
10  GW  (1  TW)  for  a  X  =  1  ym  (10  ym)  laser,  the  single  stage  energy  gain  in 
a  conventional  ICA  driven  by  a  higher-order  Gaussian  beam,  Eq.  (62),  is 
limited  to  W^  ~  5  MeV  (50  MeV).  Higher  energy  gains  require  higher  laser 
power,  and  the  effects  of  self-focusing  and  ionization  can  no  longer  be 

Q 

neglected.  Recent  experiments  at  BNL  on  the  conventional  ICA  observed  a 
3.7  MeV  energy  gain  (31  MeV/m)  of  an  injected  electron  beam  (40  MeV)  using 
a  580  MV  CO2  laser  (X  =  10.6  ym)  in  2.2  atm  of  H2  gas. 

C.  Electron  Acceleration  in  Self-Guided  ICA 

To  enhance  the  single-stage  energy  gain,  a  self-guided  ICA 
configuration  is  proposed  and  analyzed.  The  self-guided  ICA  operates  at 
laser  powers  near  the  nonlinear  focusing  power  and  intensities  near  the 
ionization  intensity.  The  self-guided  ICA  is  based  on  the  optically-guided 
higher-order  radially-polarized  optical  beam  described  in  Sec.  IV.  The 
energy  gain  is  enhanced  since  (i)  the  self-guided  ICA  operates  at  higher 
power  and  intensity  which  increases  the  accelerating  field,  and  (ii)  the 
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optical  beam  is  self-guided  which  increases  the  acceleration  length  beyond 
the  limits  of  vacuum  diffraction. 


1.  Phase  Velocity 

The  evolution  of  the  spot  size  of  the  higher-order  radially-polarized 
beam  in  the  self-guided  ICA  is  described  by  Eq.  (52).  In  addition,  the 
evolution  of  phase  velocity  v^j^  of  the  higher-order  radially- 

polarized  beam  is  given  by  Eq.  (22d)  with  p  =  1  together  with  Eqs.  (47)  and 
(51), 
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where  the  terms  on  the  right-hand  side  denote,  respectively,  the 
contributions  from  the  linear  refractive  index,  finite  spot  size,  nonlinear 
refractive  index  and  the  plasma.  For  a  self-guided  beam  with  r^  =  ’^so’ 
phase  velocity  is  <  1  and  given  by 


p  =  (l  -  0.5(Vm.„r^^)2p/pJn 


-1 

o 


(64) 


2  -1/2 

The  Lorentz  factor  associated  with  %h  -  %h'>  > 


=  ((X/ 


nn^r  )^P/Pkto  + 
o  so  NR 


1  -  n 


-2 


,-1/2 


(65) 


The  nonlinear  refractive  index  n2  can  have  a  significant  contribution  to 
the  phase  velocity. 

The  phase  velocity  can  be  controlled  by  introducing  a  small  amount  of 
background  plasma.  A  transversely  uniform  background  plasma  will  increase 
the  phase  velocity  but  have  no  effect  on  the  focusing  of  the  optical  beam. 
By  introducing  a  background  plasma,  the  right-hand  sides  of  Eqs.  (63)-(65) 
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2  2  2 

will  be  modified  by  the  addition  of  the  term  kpjj/2k  ,  where  kpj^  = 

4iiq  njj/mc  and  nj^  is  the  background  plasma  density.  The  background  plasma 
can  be  created  by  introducing  a  small  concentration  of  easily  ionized  gas, 
i.e.,  a  gas  with  a  low  ionization  energy  Uj.  In  addition,  the  background 
density  can  be  tapered  as  a  function  of  z  to  increase  the  phase  velocity 
and  optimize  electron  acceleration. 

2.  Accelerating  Gradient 

In  the  self-guided  ICA,  the  accelerated  electrons  are  acted  on  by  both 

an  axial  accelerating  E  field  as  well  as  a  transverse  field  (for  r  t 

”  z  r 

2  2 

0).  For  an  electron  near  the  axis,  r  /r  «  1,  the  field  components 

s 

associated  with  the  higher-order  radially-polarized  mode  are  given  by  E^  = 

-  (2'f2/kr  )E  sinij/  and  E  =  (■JTr/r  )E  cosV',  where  «  k^  +  0  is  the  phase, 
as  given  by  Eqs.  (44)  and  (45).  For  E^  >  0  the  electron  will  experience 
axial  acceleration  within  the  phase  regions  sintj/  >  0  and  transverse 
focusing  within  the  phase  regions  cost)/  <  0.  Simultaneous  acceleration  and 
transverse  focusing  occurs  for  ii/2  <  t  +  2iin  <  n,  where  n  =  0,±1,±2,... 

The  accelerating  gradient  E  =  (4T/ll)(X/n  r  )E  can  be  estimated  by 

zo  o  so  o 

considering  the  case  of  a  self-guided  beam  and  assuming  that  the  peak 

14  2 

intensity  is  near  the  ionization  intensity  Ij  (^10  W/cm  )  and  that  the 

power  is  near  the  nonlinear  focusing  power,  i.e.,  Ip  ^  Ij  and  P  «  For 

1/2  1/2 

this  case  we  find  that  (8JieIj/cn^)  ,  X/r^  «  (ii/2)(en^n2lj)  and 

the  accelerating  gradient  becomes 

®zo  “  <®"2^l/2n^>^^\  “  (2e/n^)(iin2li/c)^^%^^^,  (66) 
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which  is  a  function  of  only  n2  and  Ij.  For  the  equilibrium  in  Fig.  11, 
=  450  MeV. 


VI.  Ionization  and  Collisional  Losses 

In  this  section  the  attenuation  length  of  the  optical  field  due  to 
electron  collisions  and  ionization  losses  is  obtained  by  solving  the  wave 
equation,  Eq.  (1),  in  the  one-dimensional  limit  in  the  absence  of  the 
nonlinear  polarization  field  P.  Collisions  and  ionization  losses  enter  the 
wave  equation  through  the  plasma  current  density  Jp. 

In  the  presence  of  ionization  and  electron  collisions,  the  equation 
for  the  plasma  current  density  given  by  Eq.  (5)  is  modified. 


3J  /3t  =  f(o^/4nlE  -  V  J^,  (67) 

-p  V  p  m~p’ 

where  v  is  the  effective  electron  collision  frequency  for  momentum 
m 

transfer  and  is  discussed  in  Appendix  B.  Since  the  optical  frequency  is 
large  compared  to  the  collision  frequency,  Eq.  (67)  can  be 

approximated  by. 


(68) 


2 

where  Op  ~  Wp  is  given  by  Eq.  (9). 

Using  the  field  representation  in  Eq.  (6),  the  one-dimensional  wave 
equation  in  the  ^,ri  variables  becomes 


2(ik  +  3/303E/3n  =  kp  [l  -  Iv^jj/mJe, 


(69) 


where  ckp  =  Wp(^)  is  the  plasma  frequency.  Since  k  <<  |3/3^|  and 
1,  Eq.  (69)  can  be  approximated  by 

1^  E  =  -i(«Q/c)  (l  -  +  (ic/a)^)3/3^)  («pE/2),  (70) 
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where  «  =  ck,  »  =«(£)/«  is  the  normalized  plasma  frequency,  and  n  =  1 

o  pp  o  ^ 

has  been  assumed.  Substituting  E  =  E^exp(i9^)  into  Eq.  (70),  where  E^(^,h) 
and  0^(£.,n)  are  the  real  field  amplitude  and  phase,  we  obtain 


8E^/3h  =  -  r(t,1i)E^, 


(71a) 


ae 

_ c 

3h 


I  Sp(  v<=) 


aw^' 

n  Ij 

2  a^ 


(71b) 


where 

cr  =  V  -  caw^/a^  -  w$atnfE^l/9^»  ^^2) 

m  p  p  ^  P  \ 

is  the  effective  damping  rate  due  to  collisions  and  ionization.  A  similar 
result,  i.e.,  Eq.  (71),  has  been  obtained  in  Ref.  30. 

In  Eq.  (72),  the  first  term  represents  losses  due  to  collisions. 

9  .".9 

The  second  term  -  ca«^/a^  =  W«  n  /n  represents  ionization  losses  due  to 

p  p  no  p 

the  fact  that  an  electron  produced  by  ionization  in  the  presence  of  an 

optical  field  leaves  behind  a  residual  energy  approximately  equal  to  the 

oscillation  energy  e  .  The  ionization  loss  terra  can  also  be  derived  by 

OS 

2 

equating  the  rate  of  loss  of  electromagnetic  energy  -ca(E^/8Ji)/ah  to  the 

rate  at  which  electrons  are  being  produced  -caUp/a^  =  Wn^^^  multiplied  by 

2  2  2 

the  average  oscillation  energy  per  electron  =  q  E^/2m»^,  which  gives 

caE^/an  =  -  W«^(n  /n  )E^.  This  expression  for  the  ionization  losses 
0  p  no  p  o 

assumes  that  the  oscillation  energy  is  large  compared  to  the  electron 
ionization  energy,  >  Uj.  The  third  term  on  the  right-hand  side  of  Eq. 
(72)  represents  the  slippage  of  the  optical  beam  envelope  in  the  ^  =  z-ct, 
h  =  z  frame  since  the  group  velocity  v  <  c  in  the  presence  of  the  plasma. 

o 
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Three-dimension  effects  can  be  approximated  by  introducing  a  filling 
factor  f  into  Eq.  (71a).  The  rate  of  decrease  of  optical  beam  power,  due 
to  collisions  and  ionization  losses,  is  given  by 

3P/9n  = 

where 

L.  =  ck^/fv  k^  +  Wk^lf  (74) 

d  \  m  p  n) 

is  the  attenuation  length.  In  Eq.  (74),  M  is  the  tunneling  ionization 

rate,  kp  =  (4iiq^np/mc^)^^^,  k^^  =  ,  f  =  «  1  is  the 

filling  factor,  and  a  (a^)  is  the  transverse  cross-sectional  area  of  the 

generated  plasma  (laser).  As  an  example,  for  n^^  =  2.7  x  10  cm  and 
A  12  —1 

n  /n  =  10”^,  V  =  V  =  7  X  10  sec  as  discussed  in  Appendix  B. 
p  no  ’  m  en 

-2-1  3  -1 

Assuming  W/c  =  10  cm  ,  the  attenuation  distance  is  =  10  f  cm. 

The  actual  frequency  of  the  optical  beam  is  w  =  -  3(k^  +  0^)/3t  =  ck  + 

c36^/3^.  Using  Eq.  (71b)  we  find 

w/w  =  1  -  ?  3w^/3^  -  ^  3^5^/ 3^^.  (75) 

0  2  p  lo  p 

The  third  term  on  the  right-hand  side  of  Eq.  (75)  is  typically  small 

2 

compared  to  the  second  term.  Furthermore,  since  3(«)p/3^  <  0,  the  frequency 

29-31 

will  be  upshifted  as  the  optical  beam  propagates. 

VII.  Discussion 

The  propagation  of  both  fundamental  and  higher-order  Gaussian  optical 
beams  in  gases  undergoing  ionization  has  been  studied  analytically  and 
numerically.  The  propagation  model  includes  the  effects  of  the  linear  and 
nonlinear  polarization  current,  the  linear  plasma  current,  and  plasma 
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generation  via  tunneling  ionization.  Envelope  equations  describing  the 
evolution  of  the  spot  size,  amplitude,  and  phase  of  the  optical  beam  were 
derived  by  applying  the  source  dependent  expansion  method  to  the  paraxial 
wave  equation  in  the  r,  ^  =  z  -  ct/n^,  and  h  =  z  coordinate  system.  The 

envelope  equation  for  the  evolution  of  the  spot  size  r  is  given  by  Eq. 

s 

(29)  for  the  fundamental  Gaussian  beam  or  Eq.  (52)  for  the  higher-order 
radially-polarized  beam.  The  nonlinear  focusing  power  for  the  high-order 
beam  was  found  to  be  four  times  that  of  the  fundamental  Gaussian  beam, 
i.e,,  =4Pjjg  =  2X^/(Jin^n2). 

The  evolution  of  an  optical  pulse  with  an  initial  power  profile  of  the 
form  P(^)  =  P^sin^(n^/L)  for  -L  <  5  <  0  with  P^  >  Pj^  has  been  studied  by 
numerically  solving  the  envelope  equation  Eq.  (29)  (Eq.  (52)),  where  Pj^  = 
PjjG  fundamental  Gaussian  (high-order  radially-polarized) 

beam.  The  front  portion  of  the  optical  beam  where  P  <  Pj^  diffracts  up  to 
the  point  where  P  =  Pj^.  The  point  on  the  front  of  the  optical  beam  where  P 
=  P„  remains  guided,  r  =  r  .  A  narrow  region  of  the  pulse  just  behind 
the  position  where  P  >  Pj^  will  focus.  In  the  region  where  focusing  occurs 
the  intensity  and  the  ionization  rate  increase,  resulting  in  a  sharp  rise 
in  the  plasma  density.  _  The  increase  in  plasma  density  causes  the  remainder 
of  the  pulse  to  diffract.  This  behavior  is  shown  in  Figs.  3,  4,  and  12. 

Self-guided  solutions,  i.e.,  3r  /9ri  =  0,  can  result  from  a  balancing 
of  diffraction,  nonlinear  self-focusing,  and  plasma  defocusing.  These 
solutions  are  characterized  by  beam  profiles  with  P  =  P^^  at  the  front  {K  = 
0)  of  the  beam,  since  the  plasma  density  vanishes  at  the  leading  edge  of 
the  optical  beam.  As  the  plasma  density  increases  behind  the  front  of  the 
beam  (^  <  0),  the  power  in  the  self-guided  beam  also  increases  such  that  P 
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>  P„  and  3P/3^  <  0.  For  a  typical  self-guided  solution  the  amount  of 
N 

ionization  is  small,  n  /n  <  10”^.  Examples  of  self-guided  beam  profiles 

p  no  ^ 

are  shown  in  Figs.  5,  6,  and  11. 

Self-guided  optical  beams  are  subject  to  an  ionization-modulation  (IM) 
instability.  The  mechanism  of  the  IM  instability  can  be  understood  by 
considering  a  perturbation  of  the  spot  size  8r  for  which  the  power  is 
constant.  For  example,  if  the  beam  spot  size  is  increased,  i.e.,  Sr  >  0, 
the  beam  intensity  and  ionization  rate  are  reduced,  resulting  in  less 
plasma  generation.  Nonlinear  self-focusing  and  plasma  defocusing  are  no 
longer  balanced  and  the  beam  focuses.  The  focusing  beam  overshoots  its 
equilibrium  value,  i.e.,  Sr  <  0,  some  distance  behind  the  front.  When  Sr  < 
0,  the  intensity,  ionization  rate,  and  plasma  density  increase,  causing  the 
beam  to  defocus  and  again  overshoot  its  equilibrium  value.  This  focusing 
and  defocusing  of  the  beam  due  to  a  varying  amount  of  ionization  throughout 
the  beam  results  in  the  IM  instability.  For  a  fundamental  Gaussian  beam, 
the  asymptotic  linear  growth  of  the  instability  is  given  by  Sr  ~  exp[(l  ± 
i43)Ng],  where  np(0/n„o]  .  as  given  by  Eq.  (42). 

The  modulation  amplitude  increases  with  both  the  distance  from  the  front  of 
the  beam,  |^|,  and  the  propagation  distance  Vi  =  z.  The  IM  instability 
leads  to  a  disruption  (erosion)  of  the  back  of  the  beam  which  moves  toward 
the  beam  front  at  a  relative  velocity  v^  given  by  Eq.  (43).  Good  agreement 
was  obtained  between  numerical  solutions  of  the  full  envelope  equation  for 
the  IM  instability  in  the  linear  regime  and  the  analytical  expressions  for 
the  asymptotic  growth  rates,  as  shown  in  Figs.  8  and  13. 

The  results  in  Section  III  show  that  self-guiding  of  a  Gaussian  beam 
requires  a  nearly  constant  axial  power  profile  with  P  «  Pj^^  and  a  peak 
intensity  near  the  ionization  threshold.  Ip  =  Ij.  Recent  experiments  on 
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the  propagation  of  short  (200  fs,  X  =  0.8  pm)  laser  pulses  in  air  have  been 

performed  at  the  University  of  Michigan. In  these  experiments,  a  large 

portion  of  the  initial  laser  power  was  observed  to  be  confined  to  a  narrow 

spot  size  (~  40  pm)  that  propagated  some  20  m  in  air  at  1  atm.  The 

14  2 

intensity  in  the  narrow  filament  was  on  the  order  of  10  W/cm  and  partial 

ionization  was  observed  along  the  propagation  axis.  The  laser  pulse  was 

injected  into  air  with  a  large  spot  size  (~  1  cm)  and  a  peak  power  several 

times  the  nonlinear  focusing  power,  P  =  6  «  10  GW,  where  P^^^  =  1.7  GW 

for  air.  These  experimental  conditions  are  far  from  the  theoretical 

conditions  for  self-guiding  a  Gaussian  beam,  as  found  in  Section  III.  Not 

only  is  the  observed  threshold  for  nonlinear  focusing  substantially  higher 

than  the  calculated  value,  but  a  significant  portion  of  the  optical  power 

is  observed  to  reside  in  a  large  halo  surrounding  the  central  filament. 

This  large  radius  (~  1  cm)  halo  strongly  affects  the  propagation  dynamics. 

The  present  analysis  does  not  directly  apply  to  this  experiment  since  the 

observed  optical  beam  is  apparently  far  from  a  Gaussian  transverse  profile. 

In  addition,  a  self-guided  inverse  Cherenkov  accelerator  (ICA)  has 

been  proposed  and  analyzed.  In  this  acceleration  configuration  a  self- 

guided  higher-order  radially-polarized  optical  beam  propagates  through  a 

gas.  Associated  with  the  high-order  mode  is  an  axial  electric  field 

component  that  is  maximum  along  the  axis  and  can  be  used  to  accelerate  an 

injected  electron  beam.  The  phase  velocity  of  the  self-guided  mode  is  less 

than  the  speed  of  light  and  can  be  controlled  by  introducing  a  uniform 

background  plasma.  Since  the  self-guided  mode  has  a  power  near  the 

nonlinear  focusing  power,  P  =  Pj^,  the  amplitude  of  the  axial  electric 

1/2 

field  can  be  written  as  E  =  9.6(n.j/c)  I  /n  ,  where  I  is  the  peak 

ZO  Z  vl  ^ 

optical  beam  intensity.  Large  values  of  the  axial  field  amplitude. 
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typically  on  the  order  of  0.5  GV/m,  require  large  values  of  02  and  Ip. 

Since  02  is  proportional  to  the  neutral  gas  density  n^^,  the  accelerating 

field  can  be  increased  by  increasing  n^^.  To  avoid  excessive  amounts  of 

plasma  and  defocusing,  the  peak  intensity  should  be  near  the  ionization 

threshold,  I  =  In  the  self-guided  mode  the  propagation  distance  is 

P  1 

limited  by  the  IM  instability.  Since  the  number  of  e-folds  of  the  IM 
instability  scales  as  ~  Up  ,  where  the  plasma  density  np  is  a  highly 
nonlinear  function  of  Ip,  the  peak  intensity  must  be  kept  sufficiently  low. 
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Appendix  A;  Photo-Ionization  Rates 
Photo-ionization  can  take  place  in  either  the  tunneling  or  multi- 

r\  t  rt  Q 

photon  regime.  ~  These  regimes  are  characterized  by  the  Keldysh 

1/2 

parameter  Yj^  =  ,  where  U^.  is  the  ionization  energy  and  = 

(l/2)m(q  |E  |/mw)  is  the  electron  oscillation  energy.  The  Keldysh  parameter 
can  also  be  written  as  Yj^  =  where  is  the  tunneling  time,  i.e.,  the 

transit  time  of  the  electron  through  the  atomic  Coulomb  barrier.  The  low 
field  limit  (yj^  >  1)  corresponds  to  the  multi-photon  ionization  regime, 
whereas  the  high  field  limit  (yj^  <  1)  corresponds  to  the  tunneling 
ionization  regime. 

In  the  high  field  limit,  Yj^  <  1>  the  ionization  rate  can  be  determined 

by  a  tunneling  calculation  for  an  atom  in  the  presence  of  a  static  electric 

field  of  amplitude  E.  The  tunneling  ionization  rate,  i.e.,  the  probability 

25 

(per  unit  time)  of  ionization,  is  given  by 

W  =  4S^(Uj)^''^(Eg/E]exp[-  |  (Uj)^''\/e] ,  (Al) 

where  the  variables  in  Eq.  (Al)  have  been  normalized  in  terms  of  atomic 

16  —1 

constants.  In  Eq.  (Al),  =  a^c/a^  =4.1  x  10  sec  is  the 

2 

characteristic  atomic  frequency,  =  2nq  /he  =  1/137  is  the  fine  structure 

2  2  2 

constant,  a^c  is  the  characteristic  atomic  velocity,  a^  =  h  /(4e  |q |  m)  = 

5.3  X  10"^  cm  is  the  Bohr  radius,  Uj  =  Uj  is  the  ionization  energy 

in  eV,  Ujj  =  13.6  eV  is  the  ionization  energy  of  hydrogen  and  Ejj  =  Iql/a^  = 

5.2  GV/cm  is  the  hydrogenic  electric  field.  The  ionization  rate  in  Eq. 

(Al)  vanishes  at  both  small  and  large  values  of  the  electric  field,  and  has 

~  3/2 

a  maximum  at  E  =  2.3  (Uj)  Eg. 

In  the  limit  Yj^  «  1>  the  tunneling  time  is  much  less  than  the  laser 
period,  and  therefore,  Eq.  (Al)  describes  the  instantaneous 
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ionization  rate  in  the  laser  field.  For  a  circularly  polarized  laser  field 
of  the  form  E  =  |E|[cos(k:z  -  wt)ej^  +  sin(kz  -  «t)ey],  where  |e|  is 
constant,  the  ionization  rate  is  obtained  by  setting  E  =  |e|  in  Eq.  (Al). 
For  a  linearly  polarized  laser  field  of  the  form  E  =  |E|cos(k.z  -  the 

average  ionization  rate  is  obtained  by  averaging  Eq.  (Al)  over  a  laser 
oscillation  period.  The  average  ionization  rate  for  a  linearly  or  radially 
polarized  laser  field  is  found  to  be  given  by 

W  =  4(3/ll)^^2j2^(Uj)^/^Ejj/|E|)^^2exp[-  |  |E  |].  (A2) 


When  Yj^  >  ionization  occurs  by  a  multi-photon  process.  In  this 
case  photons  are  required  to  increase  the  electron  energy  by  the 
ionization  energy  Uj.  For  an  electron  ionized  from  the  ground  state  to  the 
continuum,  the  ionization  rate  can  be  approximated  by 


-2N 

W  =  A«  (2yj,)  , 


(A3) 


where  A  is  a  constant  on  the  order  of  unity  and  is  the  smallest  whole 
number  for  which  N  hft)/2n  =  U,,  i.e.,  N  is  the  number  of  photons  required 

0x0 

for  ionization. 
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Appendix  B;  Electron  Collision  Frequenc] 

As  the  electrons  in  the  weakly  ionized  gas  or  plasma  oscillate  under 
the  influence  of  the  optical  field,  they  collide  with  the  background 
electrons,  ions  and  neutral  atoms.  The  electron  collision  frequency  for 
momentum  transfer  is 


V 

m 


+  'V  f 

en’ 


(Bl) 


where  v  .  (v  )  is  the  electron-ion  (electron-neutral)  collision  frequency, 
ei  '  en' 

The  electron-electron  collision  frequency  does  not  contribute  to  because 
the  momentum  of  any  pair  of  colliding  electrons  and  associated  current 
(masses  and  charges  are  identical)  are  conserved.  Electrons  colliding  with 
electrons  will  therefore  be  accelerated  on  average  as  if  the  collision 
frequency  were  zero.  Electron-electron  collisions,  however,  lead  to 
thermalization  of  the  electrons. 

The  electron-ion  collision  frequency  is  given  by^^ 


=  4  X  lO'^l’nA^.n.Z^e  ■^''^[eV], 

€  JL  JL  OS 


where  tnA  .  is  the  Coulomb  logarithm,  (typically,  tnA  .  =  10-20),  Z  is  the 

ion  charge  state,  n^  is  the  ion  density,  n^Z  is  the  electron  plasma  density 

and  e  is  the  electron  oscillation  energy  in  units  of  eV. 
os 

33-35 

The  electron-neutral  collision  frequency  is  given  by 


<B3) 


where  a  is  the  electron-neutral  cross-section,  v  is  the  electron 
en  ® 

velocity  and  the  brackets  <  >  denote  an  averaging  over  the  electron 
velocity  distribution.  The  electron-neutral  cross-section  is  generally  a 
complicated  function  of  the  electron  velocity.  At  low  electron  velocities 
the  cross-section  is  hard-sphere-like  and  independent  of  velocity,  a 
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IS  — ^ 

where  a  -  10”  cm  is  the  hard-sphere  cross-section*  As  the  electron 
o 

velocity  approaches  the  characteristic  atomic  electron  velocity, 

polarization  scattering  is  the  dominant  process  and  a  The 

characteristic  atomic  electron  velocity  is  v^  =  a^c  where  =  1/137  is  the 

2 

fine  structure  constant  and  the  characteristic  electron  energy  mv^/2  is 

13.6  eV.  At  substantially  higher  electron  velocities  the  scattering 

becomes  Coulomb-like  and  a  ~  1/v^.  The  electron-neutral  collision 

en  e 

frequency  can  therefore  be  estimated  to  be  =  "n%'^os  '^os  ^  '^o 

v  =  n  a  V  for  v  >  v  ,  where  v  is  the  electron  oscillation  velocity, 
en  n  o  o  os  o  os 

Typically,  the  electron-neutral  collisions  is  the  dominant  collisional 

process  in  weakly  ionized  gases.  As  an  example,  consider  the  case  where 

19  —3 

the  neutral  density  is  n  =  3  x  10  cm  and  the  electron  oscillation 

no 

-3  2 

velocity  is  v  =  qE  /mo)  =  5  x  10  c  (.t  =  mv  /2  =  6.3  eV).  For  a 

OS  o  OS  os 

linearly  polarized  laser  of  wavelength  X  =  1  ym,  these  parameters 

13  2 

correspond  to  a  peak  intensity  of  I  =  3.5  x  10  W/cm  and  peak  electric 

o  _i  5  2 

field  amplitude  of  =  1.6  x  10  V/cm.  Taking  =10  cm  ,  Eq.  (B2) 

12  -1 

yields  an  electron-neutral  collision  frequency  of  =  7  x  10  sec  , 
i.e.,  an  electron-neutral  collision  time  of  ~  140  fsec. 

_3 

If  for  the  same  example,  the  gas  is  weakly  ionized  with  n^  =  10  n^^  = 

3  X  10^^  cm"^  and  the  electron  oscillation  energy  is  =  6.3  eV,  the 

electron-ion  collision  frequency  from  Eq.  (B2)  for  singularly  ionized  gas 

(Z  =  1)  is  V  .  =  8  X  10^^  sec”^,  i.e.,  an  electron-ion  collision  time  of 
'  ei 

T  .  =  =  13  psec.  For  these  parameters  the  electron-neutral  collision 

ei  ei 

frequency  is  much  greater  than  the  electron-ion  collision  frequency,  » 

V  . . 

ei 

As  another  example  consider  the  case  of  a  plasma  in  which  all  the 
atoms  are  singly  ionized,  i.e.,  Z  =  1  and  n^i^  =  n^^  is  the  initial  neutral 
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18  2 

density.  For  a  laser  intensity  of  I  =  10  W/cm  and  wavelength  X  =  1  ym, 

the  oscillation  energy  is  =  0.2  MeV.  For  a  neutral  density  of  n^^^  = 

10^^  cm”^,  the  electron-ion  collision  frequency  is  =  4  x  10^ 

sec”^  and  x  .  =  =  250  nsec. 

ei  ei 
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Appendix  C;  Radiative  and  Collisional  Losses  on  Accelerated  Electrons 

In  an  inverse  Cherenkov  or  a  laser-plasma  accelerator,  the  accelerated 
electrons  interact  with  both  the  background  electrons  and  nuclei.  This 
interaction  results  in  a  loss  of  electron  beam  energy  as  well  as  an 
increase  in  beam  emittance.  Accelerated  electrons  traveling  through 
neutral  gas  or  plasma  are  scattered  by  the  atomic  nuclei  and  emit 
bremsstrahlung.  The  emission  of  bremsstrahlung  represents  a  loss  in 
electron  energy.  In  addition,  a  highly  relativistic  electron  will  lose 
energy  due  to  collisions  with  atomic  or  plasma  electrons. 

The  rate  of  change  of  energy  of  the  accelerated  electrons  is  given  by 


dV/dz  =  , 

e  acc  rad 


-  G 


col’ 


(Cl) 


2 

where  W  =  (y  -  l)mc  is  the  electron  energy,  G  =  qE  denotes  the 

0  O-CC  2 

accelerating  gradient  and  G^.^^  ^^col^  denotes  the  radiative  (collisional) 

47 

energy  loss.  The  rate  of  electron  energy  loss  due  to  bremsstrahlung  is  ’ 


G 


rad 


(C2) 


where 

Ljj  =  ((16/3)ajn^zVfcnA]  \  (C3) 

is  the  radiation  damping  length.  In  Eq.  (C3),  n^  is  the  density  of  nuclei, 
Z  is  the  atomic  number,  =  Znq^/hc  =  1/137  is  the  fine  structure 

3.  ^ 

constant,  r  =  q^/mc^  =  2.8  x  10"^^  cm  is  the  classical  electron  radius  and 
e 

A  is  a  function  of  the  maximum  and  minimum  impact  parameters.  The  choice 

of  A  depends  on  various  factors,  such  as  the  electron  energy  range  and 

electronic  screening  effects.  For  highly  relativistic  energies,  »  W^, 

-1/3 

where  complete  screening  occurs,  A  =  233  Z^  .  At  lower  energies,  << 
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W^,  where  screening  can  be  neglected,  A  =  y.  The  critical  energy 

-1/3  2 

defining  these  regimes  is  W  =  192  Z  me  and  typically  £.nA  ~  5  -  10.  In 

C  o 

19  -3 

hydrogen  (Z  =  1)  at  1  atm  (n  =  5.4  x  10  cm  )  the  radiation  damping 

3.  3 

length  is  L|^  =  12  km,  where  we  have  taken  £.nA  =  5. 

For  a  highly  relativistic  electron,  the  energy  loss  due  to  collisions 

47  48 

with  atomic  or  plasma  electrons  is  given  by  the  Bethe  formula,  ’ 

G  T  =  4nn  r^mc^tnB,  (C4) 

coi  e  e 


where  n  =  n  Z  is  the  total  electron  density  and  B  is  the  ratio  of  the 
e  a  a 

maximum  to  minimum  impact  parameter.  The  ratio  of  radiative  energy  loss 
Eq.  (C2)  to  collisional  energy  loss  (C4)  is 


G  ,/G  T 
rad  col 


^^a  e.nA 
323  <LnB 


=  1.4  X  10 


-3 


Wg[MeV] . 


(C5) 


For  an  electron  traveling  through  hydrogen  (Z  =  1)  gas  or  plasma  the 
radiative  energy  loss  term  is  larger  than  the  collision  loss  for  energies 
above  700  MeV. 

In  addition,  accelerated  electrons  will  undergo  small  angle  scattering 

off  nuclei,  which  leads  to  an  increase  in  the  electron  beam  emittance.  The 

2  1/2 

normalized  emittance  of  the  electron  beam  is  given  by  =  yR<0  >  ,  where 

2  1/2 

R  is  the  rms  beam  radius  and  <6  >  is  the  rms  divergence  angle  of  the 
beam.  The  rate  of  growth  in  the  mean  square  angle  is  given  by 

d<eS/dz  =  16nZ^(Z^  +  Dn^r^y'^enA,  (C6) 

and  the  normalized  emittance  increases  according  to^®’^^ 

^  =  y^R^d<e^>/dz  =  R^/L  ,  (C7) 

dz  s 
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where 
5.4  X 


L”^  =  16nZ^(Z^  +  l)n^rgtnA. 
10^^  cm~^,  and  2.nA  =  5), 


For  hydrogen  at 
=  470  cm. 


atm  (Z,  =  1,  n 
a  3 
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Appendix  D;  Inverse  Cherenkov  Acceleration  vith  Bessel  (Axicon)  Beams 


In  addition  to  the  higher-order  Gaussian  optical  beam  discussed  in 

Sec.  IV,  an  inverse  Cherenkov  accelerator  (ICA)  could  be  driven  by  a  first 

7  8 

order  Bessel  optical  beam.  ’  Both  optical  beams  are  axially-symmetric, 

radially-polarized  and  have  an  axial  field  peaked  along  the  z-axis. 

Nonideal  Bessel  beams  (finite  in  transverse  extent)  can  be  formed  using 
38  39 

axicon  mirrors.  ’  Both  the  nonideal  Bessel  beam  and  the  higher-order 
Gaussian  beam  diffract,  limiting  the  acceleration  distance.  For  a  fixed 
total  optical  beam  power,  however,  the  energy  gain  in  an  ICA  can  be 
substantially  higher  when  driven  by  a  Bessel  optical  beam  as  opposed  to  a 
higher-order  Gaussian  optical  beam.  Nonlinear  self-focusing  in  the  gas  and 
the  effects  of  ionization  are  neglected  in  the  following. 

The  wave  equation  in  the  paraxial  approximation  for  a  radially- 
polarized,  axially-symmetric  field  E^.  propagating  in  a  medium  with  linear 
refractive  index  n  is 


1  ^ 
r  3r 


r"^  +  2ik 


fi_  I  i_l 

Uz  *  V  it) 


E  =  0, 
r  ’ 


(Dl) 


where  E^.  =  E^(r,z, t)e^^*^^  “^^2  +  c.c.,  E^.  is  slowly  varying  in  z  and  t 
-1  -1 

compared  with  k  and  »  ,  v  =  c/n^,  k  =  w/v,  »  =  2nc/X  is  the  frequency 

and  X  is  the  vacuum  wavelength.  An  exact  solution  to  Eq.  (Dl)  which 
maintains  a  constant  transverse  profile  is 


Ej.  =  E^Jj^(kj^r)exp  ^i(k  -  Ak)z  -  wtj 


+  c.c. 


(D2) 


where  is  the  Bessel  function  of  the  first  kind  of  order  unity,  Ak  = 

n 

k^/2k,  kj^  is  the  transverse  wavenumber  and  E^  is  the  radial  field 

amplitude.  The  ideal  Bessel  field  in  Eq.  (D2)  (infinite  in  transverse 

50 

extent)  is  nondiffracting  in  the  sense  that  the  transverse  profile 
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remains  constant.  The  power,  however,  contained  within  an  ideal  Bessel 

2 

beam  is  infinite  since  JdrJ^(kj^r)r  =  ®  when  integrated  from  r  =  0  to  r  =  ®. 
Associated  with  the  radially-polarized  field  in  Eq.  (D2)  is  the  axial  field 


ik,  E 

F  -  ^  ° 

”  (k  -  Ak) 


J^(kj^r)exp  ^i(k  -  Z\k)z  -  «tj 


+  c.c. 


(D3) 


The  axial  accelerating  field  in  Eq.  (D3)  is  peaked  along  the  z-axis  and  has 
axial  phase  velocity  v^j^  =  «/(k  -  Ak), 

Vph  v[l  +  (kj^v/w)^/2] ,  (D4) 

2 

which  can  be  less  than  c  for  n^  >  1  +  (kj^c/co)  /2,  where  it  is  assumed  that 

(kj^c/«)^  «  1. 

The  ideal  Bessel  beam  consists  of  an  infinite  number  of  rings  (lobes) 

extending  radially  to  infinity  and  having  a  radial  width  of  r^  ^  n/k^* 

Since  the  asymptotic  form  (kj^r  »  1)  for  the  Bessel  function  is  J^(k^r)  ^ 

(2/jtkj^r)^'^^cos(kj^r  -  3ii/4),  the  power  in  each  ring  is  essentially  the  same. 

If  the  power  in  each  ring  is  denoted  by  P^,  then  the  total  power  contained 

in  a  nonideal  Bessel  beam  of  a  finite  radial  extent  is  P  NP^,  where 

N  =  R  /r  is  the  number  of  rings.  In  principle,  the  number  of  rings  can 
max  o 

be  large,  N  »  1.  A  nonideal  Bessel  beam  consisting  of  N  rings  diffracts 

51 

away  sequentially  starting  with  the  outermost  ring.  The  outermost  ring 

2 

diffracts  after  a  distance  ~  the  next  ring  diffracts  after  a 

distance  2iir^/X,  and  so  on  until  the  innermost  ring  diffracts  away  after  a 

distance  ~  Nnr^/X.  Hence,  the  maximum  propagation  distance  of  a  nonideal 
o 

Bessel  beam  consisting  of  N  rings  of  width  r^  is 


L 


max 


=  NZ 


Ro’ 


(D5) 
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where  Z_  =  Jir^/X  is  the  Rayleigh  length  associated  with  the  individual 

KO  O 

rings,  assuming  n^  =  1. 

The  maximum  energy  gain  in  the  ICA  driven  by  a  nonideal  Bessel  beam  is 


W  =  -qE  L  , 
max  ^  zo  max’ 


(D6) 


assuming  that  the  axial  phase  velocity  is  matched  to  the  electron  velocity, 

where  E  =  (k,/k)E  is  the  axial  accelerating  field  along  the  z-axis  given 
zo  1  o 

by  Eq.  (D3).  The  radial  field  amplitude  in  terms  of  the  power  within  a 
ring  is 

1  /9 

E  =  (2n/r  )(P^/cn^)  .  (D7) 

Using  Eqs.  (D5)  and  (D7),  the  maximum  energy  gain  from  Eq.  (D6),  in  terms 

1/2  1/2  2  1/2 

of  the  total  optical  power,  is  W  =  C  N  P  ,  where  C  =  -qit  /(cn  ) 

In  practical  units, 

W  (MeV)  =  1.7  N^^^[P(GW)]^^^,  (D8) 

m^x 


for  an  ICA  driven  by  a  nonideal  Bessel  beam.  If  a  higher-order  Gaussian 
optical  beam  of  the  same  total  power  P  were  used  instead  of  the  nonideal 
Bessel  beam,  the  maximum  energy  gain,  Eq.  (62),  would  be  “ 

I/O  1/2 

2.3[P(GW)]  .  The  energy  gain  in  the  ICA  is  ~  N  times  greater  for  a 

nonideal  Bessel  beam  as  compared  to  a  higher-order  Gaussian  beam  of  the 
same  total  power. 

The  ratio  of  the  accelerating  gradient  for  the  nonideal  Bessel  beam, 

^zB’  higher-order  Gaussian  beam,  E^^,  is  E^g/E^g  ~ 

/r  )^,  where  r  is  the  spot  size  for  the  higher-order  Gaussian 
'  so  o  so 

beam. 
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Fig.  1  —  Schematic  of  a  fundamental  Gaussian  optical  beam  propagating  in  a  gas  undergoing  ionization.  When  P  >  Pnq,  the  beam  self¬ 
focuses  and  the  intensity  increases  causing  ionization.  Plasma  is  generated  along  the  beam  axis,  tending  to  defocus  the  beam.  Self-guiding 
results  by  balancing  nonlinear  self-focusing  with  plasma  defocusing. 
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Fig.  2  —  Surface  plots  of  (a)  optical  pulse  intensity  I  and  (b)  plasma  density  np  plotted  versus  (r,^)  at  j;  =  z  =  0  for 
a  fundamental  Gaussian  beam  propagating  in  air  at  1  atm,  propagation  is  towards  the  right.  Plot  (c)  shows  the 
intensity  I  (solid  line)  and  plasma  density  np  (dashed  line)  versus  r  at  the  pulse  center  (X  —  -30  jtim). 


Fig.  2  (continued)  —  Surface  plots  of  (a)  optical  pulse  intensity  I  and  (b)  plasma  density  rip  plotted  versus  (r,$)  at 
=  z  =  0  for  a  fundamental  Gaussian  beam  propagating  in  air  at  1  atm,  propagation  is  towards  the  right.  Plot  (c) 
shows  the  intensity  I  (solid  line)  and  plasma  density  np  (dashed  line)  versus  r  at  the  pulse  center  (^  =  —30  ^m). 
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Fig.  2  (continued)  —  Surface  plots  of  (a)  optical  pulse  intensity  I  and  (b)  plasma  density  np  plotted  versus  (r,0  at 
^  =  z  =  0  for  a  fundamental  Gaussian  beam  propagating  in  air  at  1  atm,  propagation  is  towards  the  right.  Plot  (c) 
shows  the  intensity  I  (solid  line)  and  plasma  density  np  (dashed  line)  versus  r  at  the  pulse  center  ($  =  -30  /xm). 
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Fig.  3  (continued)  —  Spot  size  Fj  (dashed  line),  intensity  I  (solid  line),  and  plasma  density  np  (dotted  line)  plotted 
versus  ^  on  axis  at  (a)  =  z  =  0,  (b)  6  cm,  (c)  8  cm,  and  (d)  10  cm.  The  initial  peak  intensity  is  Ipo  =  3.0  x  1 

W/cm^,  the  initial  spot  size  is  t^q  =  80  ^m,  and  the  peak  power  is  Pq  =  3.0  GW  =  1.1  Png-  The  direction  of 
propagation  is  towards  the  right. 


Fig.  6  —  Equilibrium  profiles  of  power  P  (solid  lines)  and  plasma  density  np  on  axis  (dashed  lines  plotted  versus  ^ 
for  three  different  values  of  leading-edge  (^  =  0)  intensity  Ip  =  Ij  =  5.0  x  10^^  W/cm^,  and  I2  =  5.1  x  10^^ 
W/cm^,  and  I3  =  5.2  x  10*^  W/cm^. 
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Fig.  7  (continued)  —  Spot  size  Fj  (solid  line)  and  plasma  density  np  on  axis  (dashed  line)  plotted  versus  ^  at 
(a)  =  z  =  0,  (b)  400  cm,  (c)  450  cm,  (d)  500  cm,  (e)  550  cm,  and  (0  600  cm,  for  an  initially-matched  laser 

pulse  with  Ip  =  3.0  x  10*^  W/cm^  and  Fjd)  =  78  fim  propagating  in  air.  The  direction  of  propagation  is  towards 
the  right. 


Fig.  7  (continued)  —  Spot  size  (solid  line)  and  plasma  density  np  on  axis  (dashed  line)  plotted  versus  ^  at 
(a)  7;  =  z  =  0,  (b)  400  cm,  (c)  450  cm,  (d)  500  cm,  (e)  550  cm,  and  (0  600  cm,  for  an  initially-matched  laser 
pulse  with  Ip  =  3.0  X  10^^  W/cm^  and  78  (im  propagating  in  air.  The  direction  of  propagation  is  towards 


Fig.  8  —  Perturbed  radius  log(|  Ar|)  (solid  line)  and  number  of  e-folds  (dashed  line)  plotted  (a)  versus  $  at  fixed 
77  =  550  cm  and  (b)  versus  tj  at  fixed  ^  =  -49  ^tm.  Here,  Ar  =  -  T^Q)/r^  is  determined  from  the  integration  of 

the  envelope  equation,  Eq.  (29),  while  Ng  is  given  by  Eq.  (42). 
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Fig.  8  (continued)  —  Perturbed  radius  log(|  Ar|)  (solid  line)  and  number  of  e-folds  Ng  (dashed  line)  plotted  (a)  versus 
^  at  fixed  7j  =  550  cm  and  (b)  versus  ij  at  fixed  $  =  -49  ixm.  Here,  Ar  =  (r^  -  rso)/rso  is  determined  from  the 
integration  of  the  envelope  equation,  Eq.  (29),  while  Ng  is  given  by  Eq.  (42). 
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Fig.  9  (continued)  —  Surface  plots  of  (a)  optical  pulse  intensity  I,  (b)  plasma  density  np,  and  (c)  accelerating  field 
=  I  E^l  plotted  versus  (r,0  at  tj  =  z  =  0  for  a  higher-order  radially-polarized  beam  propagating  in  H2  at 
30  atm.  The  direction  of  propagation  is  towards  the  right.  Plot  (d)  shows  the  intensity  (solid  line),  accelerating  field 
Ej  (dashed  line),  and  plasma  density  rip  (dotted  line)  plotted  versus  r  at  the  pulse  center. 
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Fig.  9  (continued)  —  Surface  plots  of  (a)  optical  pulse  intensity  I,  (b)  plasma  density  tip,  and  (c)  accelerating  field 
=  lE^I  plotted  versus  (r,|)  at  17  =  z  =  0  for  a  higher-order  radially-polarized  beam  propagating  in  H2  at 
30  atm.  The  direction  of  propagation  is  towards  the  right.  Plot  (d)  shows  the  intensity  (solid  line),  accelerating  field 
§3  (dashed  line),  and  plasma  density  np  (dotted  line)  plotted  versus  r  at  the  pulse  center. 
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Fig.  9  (continued)  —  Surface  plots  of  (a)  optical  pulse  intensity  I,  (b)  plasma  density  rip,  and  (c)  accelerating  field 
£  =  |E^|  plotted  versus  (r,^)  at  rj  =  z  =  0  for  a  higher-order  radially-polarized  beam  propagating  in  H2  at 
30  atm.  The  direction  of  propagation  is  towards  the  right.  Plot  (d)  shows  the  intensity  (solid  line),  accelerating  field 
(dashed  line),  and  plasma  density  np  (dotted  line)  plotted  versus  r  at  the  pulse  center. 


r  =  x^J^|2  (dotted  line)  plotted  versus  ^  at  (a)  ij  =  z  =  0,  (b)  60  cm,  (c)  75  cm,  and  (d)  90  cm,  for  an  initially 
matched  optical  pulse  with  Ip  =  3.2  x  10*^  W/cm^  propagating  in  H2  at  30  atm.  This  initial  radius  is  r^C^  =  37 
The  direction  of  propagation  is  towards  the  right. 


Fig.  12  (continued)  —  Accelerating  field  on  axis  (solid  line),  spot  size  (dashed  line),  and  plasma  density  np  at 
r  =  (dotted  line)  plotted  versus  ^  at  (a)  tj  =  z  =  0,  (b)  60  cm,  (c)  75  cm,  and  (d)  90  cm,  for  an  initially- 

matched  optical  pulse  with  Ip  =  3.2  x  10*^  W/cm^  propagating  in  Hj  at  30  atm.  This  initial  radius  is  rs(0  =  37  ^ 
The  direction  of  propagation  is  towards  the  right. 


Fig.  13  —  Perturbed  radius  log(|Ar|)  (solid  line)  and  number  of  e-folds  Ng  (dashed  line)  plotted  (a)  versus  ^ 
at  fixed  17  =  z  =  75  cm  and  (b)  versus  z  at  fixed  ^  =  —40  fim.  Here,  Ar  =  (r^  —  rso)/rso  is  determined  from 
the  time  integration  of  the  envelope  equation,  Eq.  (52),  while  is  given  by  Eq.  (58) 
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Fig.  13  (continued)  —  Perturbed  radius  log(|  Ar|)  (solid  line)  and  number  of  e-folds  Ng  (dashed  line)  plotted 
(a)  versus  ^  at  fixed  t/  =  z  =  75  cm  and  (b)  versus  z  at  fixed  ^  =  -40  fim.  Here,  Ar  =  (r^  -  r,o)/rso  is 
determined  from  the  time  integration  of  the  envelope  equation,  Eq.  (52),  while  is  given  by  Eq.  (58). 


